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FOREWORD 


This report presents the fina! results of one of the 46 projects comprising the military-effects 
program of Operation Plumbbob, which inclured 24 test detonations at the Nevada Test Site in 
1967. 

For overall Plumbbobd military-effects information, the reader is referred to the “Summary 
Report of the Director, DOD Test Group (Programs 1-9),” ITR—1445, which incindes: (1) a 
deseription of each detonation, including yield, zero-point location and environment, type of . 
device, ambient atmospheric conditions, etc. ; (2) a discussion of project results; (3) 2 summary 
of the objectives and results of each project; and (4) a listing of project reports for the military- - 
effects program. 


\ ABSTRACT 


Project 1.4 objectives were to measure underground effects of 2 nuciear air burst (Shot Pris- 
cilia; 36.6 kt) as they vary with time, depth, and ground range, particularly in the region of 
high pressure; furnish, from these measurements, input data to other projects; and analyze 
there measurements with results of other testy. At 750 and 1,050 [eet from ground zero, 
acceleration, stress, and strain measurements were made at several depths down through 50 
fect, including two measurements cach of horizontal acceleration and stress. At 450, 550, 
650, 750, 850, 1,050, and 1,350 feet, vertical acceleration and stress were measured at 50 
and 10-foot depths. 

Recards were obtalged oa 52 out of 64 gage channels installed. Losses were caused by 
mechanical failure of one oscillograph and by misceHaneous individual channel failures. 

Wave forms of acceleration and velocity showed no ideal or classical shape but could be 
grouped in six categories according to their characteristics. 

Qutrunning occurred at the ground surface at 2,500-foot ground range from a signal origi- 
nating at 1,900-foot ground range. However, outrunning can occur at closer ranges for deep 
measurements, and refracted signals may ce recorded after arrival of local effects, as evi- 
denced by acceleration measurements. 

Attemzation of maximum downward acceleration at 5- and 10-foct depths varied between 30 
and 45 percent except at 550- and 650-foot ranges where it was negligible. At greater depths, 
wave theory concerning energy transfer at an interface between two materials was borne out. 
Horizontal (outward) acceleration at 10- and 50-foot depths was attenuated less with depth than 
was the corresponding peak downward acceleratior. 

Peak downward velocity followed an exponential decay law rather than a power law decay 
characteristic of downward' acceleration. At 275~psi level, horizontal (outward) velocity showed 
somewhat less attenuation vith depth than the downward component. At the 100-psi level, peak 
outward velocity at £0-foot depth was twice that at 10-foot depth, owing to signals from sources 
closer to ground zero. 

Attemation of peak displacement corresponded closely to attenuation of peak velocity. 

Attenuation of maximum vertical stress was slight between the surface and 5-foot depth, and 
stress decreased by half for every 10-fooct increase in depth, except at 50-foot dcpth where it 
increased. Stress measurements on this project were not considered entirely successful, de- 
spite extreme caution exercised in gage placement and backfill procedure. 

At 275-psi overpressure, peak strain decreased abruptly between 1- and 20-foot depths, 
leveling off to 2 constant value at greater depths. At 100-psi overpressure, vertical strain 
showed almost no change with depth This difference between the two stations could probably 
be traced to the longer rise time of the overpressure at 100 psi than at 275 pst. 

The velocity-jump peak overpressure ratio increased with decreasing pressure, with no 
apparent systematic variation with yleld, overpressure level or wave form, or test area. 
Experimental ratios agreed well with the theoretical result. Peak vertical displacement- 
overpressure tmpulse ratio data were too scattered to allow firm conclusions. 

From displacement-rcsponse spectra, the change in the character of the response appeared 
to be associated with the interference of the refracted ground-transmitted wave and not with the 
local ground wave. Normalized velocity spectra for 5-, 10-, and 50-foot depths showed similar 
maxima although the {requency at which this maximum occurred decreases with increasing depth. 

In stress-strain relations, it was tentatively concluded that laboratory triaxial tests were 
more useful in correlation with blast results than were compaction tesis. 
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Chopter | 
INTRODUCTION 


1.1 OBJECTIVES 


The objectives of Project 1.4 were to: (1) measure underground effects of a nuclear air burst 
as they varied with time, depth, and ground range, particularly in the region of high pressure; 
(2) use these measurements to furnish input data to other projects; and (3) analyze these data in 
combination with the results of other tests to establish criteria for the prediction of underground 
effects. 

The quantities measured were earth acceleration, earth stress, and earth strain. 

In addition, to obtain more information on the test medium, seismic studies were performed 
at Frenchman Flat. 


1.2 BACKGROUND 


Underground effects of a nuclear air Furst, as opposed to underground effects of an uncer- 
ground burst, have not been extensively investigated. The most compiete full-scale study of 
these effects was on Operation Upshot-Knothole (Reference 1), when detailed measurements of 
earth acceleration, stress, and strain were made at one depth and at one ground range on two 
shots and when measurements of earth stress were made at three depths at several ground 
ranges on the same shots. On Operation Tumbler (Reference 2), earth acceleration was meas- 
ured on four shots; measurements were made at two depths on one shot and at three depths at 
‘two ground ranges on all four shots. Only one of these shots was in Frenchman Fiat; the other 
three were in the Yucca Flat T7 area, where geological conditions were considerably different. 
All these measurements were made at pressure levels considerably lower than the regions of 
present interest, Small-scale studies of similar phenomena conducted by Stanford Research 
Institute (SRD on Project Mole (Reference 3) and by SRI on the high-explosives series of Opera- 
tion Buster-Jangle (Reference 4) provided little useful information; the limited frequency re- 
sponse of the instruments used resulted in under-registration of the phenomena. Therefore, 
extrapolation of these small-scale results to full-scale effects was prohably not useful. 

Rather extensive small-scale and full-scale studies of underground effects of underground 
high explosives have been conducted, but the transmission and loading mechanisms under high- 
explosive test conditions differed so markedly from those of the conditions of this project that 
extension of their results to these problems was difficult if not impossible. 

Requirements have been set up for criterla for protective construction against large air and 
surface bursts in the pressure regions of 100 psl and above. Such requirements imply that, 
wherever possible, underground construction must be used for maximum protection. It is be- 
lieved that most of the loading of underground structures in these pressure regions from such 
shots will be produced by the air-biast slap on the surface in the immediate vicinity of the struc- 
ture (local effects), rather thin by energy transmitted through the ground from the regions 
closer to ground zero (remote effects). This project was established to obtain quantitative data 
on underground phenomena in these regions and to explore these phenomena ‘sufficiently to per- 
mit at least a tentative application of the results to similar effects in soils other than those 
characteristic of Frenchman Flat. Because of limitations {mposced at the Nevada Test Site (NTS), 
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tt was not possible to plan a large surface shot; therefore these studics were based on the effects 
of a moderately large (35 to 40 kt) ater burst. 

A compicte description of the ground motion Induced by air blast Lrom 2 high-explosive or 
nuclear detonation requires spatial and temporal specification ef toth the horizontal radial 
(hereafter referred to simply as the horizontal componcnt) and vertical components of: (1) 
particle acccleration; (2) particle velocity; (3) particle dispiacement; and (4) earth strain. 
Theoretically, only one of these need be presented alone with earth stress since any onc of the 
first three may be derived from the other by differentiation or integration with respect to time, 
and strain may be found by differentiation of displacement with distance and vice versa. To- 
gether these parameters form a redundant set of data since the relationship between them is 
simply one of geometry. On the other hand, if the medium exhibits Inclastic behavior, the 
relationship of stress to particic velocity, for cxample, requires a great deal more knowledge 
than {fs now available. Furthermore, if the medium is strain-rate sensitive (as soll is suspected 
to be) or passesses visco-elastic properties, then stress depends not only on the instantancous 
values of strain and/or particle velocity but on their past history as welt. 

From a practical standpoint, specification of ground motion 2s a single parameter [s seriously 
limited by the small number and low density of observations made on any one experiment and also 
by the inherent complexity of soil-particle motion compared with that of air. At present, ground 
motion must be specified as positive and negative peak values, time durations, and representa- 
tive pulse shapes. Hence, salient features of acceleration, velocity, displacement, and strain 
may be lost. If these salient features are required, as in computation of response spectrum 
(see Chapter 4), it is necessary to refer to the original records. 

In Seeking 2 valid correlation of ground motion data with air blast input data, answers to 
Several fundamental questions are required. These may be set forth as follows: 

1. What relationships are to be expected between the various elements of ground motion and 
air-blast input? How are these relationships modiiied by soil properties and by variation of 
Soil properties with depth? 

2 Are the abcve relationships expected to vary with yield of the device, ground range, 
and/or stress level? If so, are these variations of major importance? 

3. For a particular test area, Le., a relatively fixed set of soil properties, are data from 
various experiments internally consistent? 

4. Do systematic differences in the correlations appear when data from one test area are 
compared with those from another? Are these differences consistent with known variations af 
soil properties? : 

Obviously the answers to Items 1, 2, and 3 depend a great deal on theoretical knowledge of 
the problem since the number of measurements involved are not sufficient to form an independ- 
ent empirical evaluation of all factors involved. 


1.3 THEORY 


Adaption of historical theories of the propagation of stress through soil is limited in that 
either (1) flows are slow under applied loads or stresses (soil statics), or (2) stress waves are 
of low intensity (exploration and carthquake seismology). The present need is for thcorics and 
Special experiments that consider stress waves of high intensity and rapid rates of loading. 


1.3.1 Models for Soil Reacting to Stress. In this section are reviewed the well-known math- 
ematical models (Figure 1.1) for solids which are subjected to stress, and the limitations and 
range of applicability of these models are described. 

Linear Elasticity: Statics. For the study of small strains, linear elasticity may 
be a satisfactory model (Column 1 of Figure 1.1). In this model, soil is assumed to act essen- 
tially as a linear spring supporting a mass. Because of the success of this model in studying 
complicated problems involving steel, it has received a tremcndous 2mount of attention. 

The theory of elasticity without further qualification may be understood to mean 2 theory 
based on the assumption that there is a unique relationship between stress and strain; more 
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explicitly, that for the material under consideration, cach component of stress can be caicu- 
lated solely from 2 knowledge of the local strain. This assumption is believed to be false for 
all known materials; but it is a useful approximation. In actual practice, the previous history 
of the material, the rate of strain, and even the derivative of the rate of strain can aficct the 
values of the stress at cach point. Spectacular examples of dependence on these quantities arc 
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Figure 1.1 Typical stress and strain diagrams for various 
models of soil behavior (sudden stress). 


given by bouncing putty, bread dough, natural rubber, and synthetic plastics. 

The complications or modifications of pure (linear) elastic theory which must be considered 
can be separated into two classes (Table 1.1). The first class is the set of complications in 
space, e.g., vertical anisotropy. The second class is the set of complications which are in- 
trinsic in every cubic millimeter of the sotl. Certainly the two types uf complications can be 
combined: viscosity of soil can change with depth, etc. 

Nonlinear Elasticity. It is possible to conceive of a substance in which the Lame’ 
constants, or (since they are related) in-which Young’s modulus and Poisson’s ratio, depend on 
the strain. A simpic example of such a substance is a spring which becomes stiffer in compres- 
sion. A helical spring will exhibit this characteristic when it closes coil to coil. 

For 2 spring which becomes stiffer with increasing strain, the mathematical model is that of 
a substance which consists of various fibers, some of which offer no resistance to strain until 
the strain reaches a finite value. A corresponding mode} is possible for a spring which becomes 
less stiff with increasing strain. These models can even be combined to give a model of a sub- 
stance, the stiffiness of which first decreases and then increases with increasing Strain. For all 
these models, the supposition is that the strain is completely recovered when stress is relieved. 
These models can probavly be used successfully for certain analyses; since they are tou simple 
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te explain certain other experimental results, they scem not to have been used cxtensively In the 
literature, 

Plasticity. By nonlinear behavior is meant usually 2 medium which has properties diff- 
erent even from those of the paragraph above. If a solid undergoes permanent deformation or 
set when it is strained, it is said to be plastic (clay, for example). A schematic stress-strain 
relation for a plastic material is that of Figure 1.2. 

This diagram is similar in some respects to Figure 1.3 which is a stress yorsus strain curve 
determined from dynamic measurement: obtained In typical silty clay soil. 

In Figure 1.2, the matcrial is called pericctly plastic if AB {ts 2 horizontal line and OA and 
BC are paraliel straight lines. Such a material behaves as a linear elastic material with no per- 


TABLE L1 REFINEMENTS OF LINEAR ELASTIC THEORY 


Compllcationa 


Sopeliations inc Intrinsic to the Medium 


Space-Dependert Velocities Nonltnearity 
Interfaces Viscosity 

Free Surfaces Visco-Elasticity 
Layering Intrinsic Anisotropy 
Refractions {as in crystals) 
Reflections 

Interference 

Lenses 

Bealders 

Inhomogencities 


manent set if the stress remains below a yield value; the stress can never exceed this yield 
value; when the yield value is reached, the material flows plastically. A more general curve, 
such as that of Figure 1.2 would be obtained by postulating a model with several fibers; ezch 
fiber is perfectly plastic, but the fibers have different yield values. 

Viscosity. It is convenient to study ihe propagation of waves in a plastic suostance by 
assuming that 2 certain amount of strain energy is dissipated at the wave front. One possible 
mode of dissipation, which seems physically reasonable, is that the rate of loss of energy de- 


STRESS 


Figure 1.2 Schematic stress-strain relation 
for 2 plastic material. 


pends on particle velocity; this is viscous friction. The mechanical analog for this mode is that 
of 2 spring and dashpot in parallel (see Column 2 of Figure 1.1). For this model, the attenuation 
per cycle is a function of the frequency. 

Visco-Elasticity. If stress is applied to a sample of soil rapidly and then suddenly 
removed, the graph of strain against time is like that in Columa 3 of Figure 1.1. A plausible 
model is that at the top of the column. The upper spring and dashpot allow for permanent set 
and rapid relief of stress; the spring and dashpot in parallel allow for gradual relief of the re- 
maining portion of stress. Indeed, any number of parallel spring-dashpot units (Voigt units) 
can be connected in series to make 2 still more complicated model. The model in question can 
be thought of as a system of three Voigt units in series. The dashpot constant in the first unit 
is zero: it is effcetively a spring alone. The spring constant in the second unit is zero: it is 
effectively « dashpot alone. In the third unit, the spring constant and the dashpot constant are 
both different from zcro. 
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Solld Friction. Ina solid friction modcl, (Column 4, Figure 1.1) dissipation of energy 
is assumed to depem! on particle displacement but not on particle velocity. The attenuation per 
cycle Ls thus independent of the frequency. Certata laboratory experiments of Born (Reference 
5) on oven-dried rocks support this modcl. For wet rocks, Born shows that viscous as well a3 
solid friction damping is present. This result will probably be true for clay soll as well. In 
the model for solid friction used by Born, the damping [orce Is taken to be proportional to dis- 
placement from the nevtral position, This is not stated by him, but can be inferred from the 
formulas which he uses. This kind of damping is also termed structural damping. 

Other Mathematical Models. Besides the possibility, already alluded to, of using 
more than three Voigt units in series, the above schemes are not the only ones that have to be 
considered in making mechanical] models of the behavior of soil. Unfortunately soil is not iso- 
tropic. The most common variation in properties which it is necessary to take into account is 
horizontal layering. For example, there can be horizontal layers, in whlch seismic velocities 
are widely different. It is frequently truc also that in any one stratum the horizontal velocity 
is not the same as the vertical velocity. 

The propagation of waves in a layered medium is 2 special topic. Some wave paths which 
enter the medium can be refracted and returned to the surface. Even the most straightforward 
problems in wave propagation in a layered (linear) elastic medium are difficult. 

Table 2 gives practical limitations of some of the theories which are disccssed here. The 
table shows two thirgs: (1) no single, simple model will explain all the phenomena observed 
when shock waves pass through soll; and (2) the model to use in studying particular types of 
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Figure 1.3 Actual record of stress versus 
Strain in typical silty clay. 


problems must be a compromise between the faithfulness with which one wants to portray all 
elements of soil behavior and the economics of producing many numerical solutions. This is to 
say that the least complicated model which will give reasonably good agreement with experimen- 
tal facts is the model to use. However, some rudimentary numerical solutions will always be 
necessary before one can choose the simplest and yet the most appronriate model for a given 
problem. 
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1.3.2 Homogencous Elastic Solid. Within the framework of present theoretical knowledge, 
several postulated relationships derived from analysis of a homogencous clastic solid can be sct 
dawn (Reference 6). 


The assumptions under which this analysis was derived are as follows. First, the plan radius 
of curvature of the air-blast wave is large compared with the depth under investigation so that in- 
duced ground motion may be regarded as two-dimensionaL Second, it is assumed that the air- 


TABLE 12 PHENOMENA EXPLAINED BY VARIOUS MODELS OF SOIL BEHAVIOR 


Dissipation Permancnt Pee a Dispersion of 
Model of Energy Set Wave Paths Surface Waves 
to Surface 
Linear Elastic (Isotropic medium) No No No No 
Nonlinear Elastic 
@sctropic medium) No No Yes* Yes* 
Elastic with Vertical Anisotropy, 
(ncluding discontinuous layers) Xo No Yes Yes 
Solid (structural) Friction Yes - No No No 
Viscosity Yes No No No 
Visco-Elastic Yes Yes No No 


* These entries imply possible explanation of the phenomena, especiaily in connection 
with waves of high stress. 


blast wave of invariant magnitude is moving at a constant velocity, U, and that sufficient time 
has elapsed so that a steady ground disturbance pattern appears to an observer fixed im the coor- 
dinates of the moving blast wave. 

An elastic solid exhibits two wave speeds: one corresponding to the propagation of waves af 
dilatation (or compression) Cz, the seismic L (longitudinal) wave velocity, and another corre- 
sponding to propagation of components of rotation, Cg, the seismic S (transverse, shear) wave 
velocity. Three distinct theoretical cases are present in the above situation, 


1. U<Csg< Cy, the subseismic case. 
2 Cg< U< Cy, the transseismic case. 
3. Cg< Cy < U, the superseismic case. 


The theory in Cases 1 and 2 says that the ground motion outruns the air-blast wave; the earth 
moves downward and away from ground zero. Although ground motion outrunning is observed in 
the field, the initial motion of the earth is generally upward and away from ground zero. This 
motion is the result of refraction of the wave front due to the vertical gradient of seismic veloc- 
ity. The overpressure at which outrunning occurs depends on the device yield and the seismic 
velocity gradient (Section 1.3.3). For the yields and seismic velocity gradients at the NTS (French- 
man Fiat), this overpressure ts approximately 10 psi. 

At higher overpressures, the blast-wave velocity is large enough so that it moves superseis- 
mically (Case 3) and the ground motion lags progressively behind that at the surface as in Figure 
1.4. At the surface, the vertical ground motion and the blast wave are related simply, viz., 


1 aP ij 
v= sa- HR * *(&) 


P U 
wv = sar °° (ez) 
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Where Ay, Uy, and Dy arc the vertical particle acceleration, vertical particle velocity, and 
vertical particie displacement; p is the soll density; P ts overpressure; and t istime. The 
function, + (derived from results in Reference 6), is shown in Figure 1.5 and is nearly unity 

for U/Cy, > 1.5. At small depths, one would expect the peak negative displacement to be pro- 
portional to the positive overpressure impulse and the peak negative velocity to be proportional 
to the peak overpressure, Py. (The usual convention is preserved; positive motion is upward 
and outward faway from ground zero] and negative motion is downward and inward [toward ground 
zeroj.) The peak acceleration depends on the overpressure wave form and for an ideal shock 
wave would be infinite. In reality, all shocks have a finite rise time, tee so chat 


Ay = 


pCy ss tr 


If one postulates that all nonprecursor shocks have approximately the same rise time, then 
one might expect the acceleration to be also proportional to the peak overpressure. 

These simple statements form the basis of correlation of near-surface ground motion. At- 
tenuation of ground motion with depth must at present be treated empirically. 


1.3.3 Effects from Remote Sources. When the propagation velocity of the air blast is larger 
than the compressional (L) wave velocity of the ground, signals will not be observed underground 
prior to arrival of the wave front (Figure 1.4). In a homogeneous medium, the wave front will 
be bent upward at increased depths (Curve a, Figure 1.4); however, if the seismic velocity in- 
creases with depth, the wave front will tend to be bent downward (Curve b, Figure 1.4). 

As the air-blast velocity decreases with distance from ground zero, the information given 
the earth will eventually arrive prior to arrival of the direct signaL. This may happen either 


Oca mens 


Figure 1.4 Schematic diagram of superseismic wave front. 


when 2 near-surface seismic velocity is greater than the air-blast velocity or by ground trans- 
mission along a curved path which dips into the higher velocity, lower earth strata, Le., by 
seismic refraction. The former corresponds to the transseismic case (Case 2), and the initial 
signal is down amd outward as in the superseismic case. The cnset of acceleration is, however, 
more gradual. The latter is distinguished by an upward and outward first signal and Is by far 
the predominant type of outrunning observed at both the NTS and the EPG. 

The outrunning signal arriving at a gage by refraction is generally called the air-induced, 
ground-transmitted wave, or, simply, the remote effect. The signal arriving by the direct 
path as the air blast passes over the yage is called the air-induced, direct wave, or, simply, 
the local effect. It is a misnomer to describe only the latter as the 2ir-blast-induced wave 
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since both classes of motion are air-blast induced, with the only difference being the point of 
origin and the transmission path. (Since the air-induced remote effect ly duc to an extended 
source, It should not be confused with the dircctly-induccd ground motion resulting from an 
undergrounc detonation.) 

Qutrunning of the ground wave has considerable importance in that the character cf the ground 
motion chioges when this occurs. Ifa scismic refraction survey has becn made of the area, the 
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Figure 1.5 versus U/C, and Poisson’s ratio. 


technique of constructing the ground motion arrival time distance is simple (Figure 1.6). 

In the lower portion of Figure 1.6 is shown the measured refraction survey of Frenchman 
Flat made during 1957. The concave portion of the refraction arrival time-curve is due to the 
increase of seismic velocity with depth. The air-blast arrival curve is first constructed. At 
each point along this curve 2 signa! is generated which travels outward at the rate prescribed 
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Figure 1.6 Construction of ground motion arrival 
time curves, Tumbler Shot 1, 5-foot depth. 


octane 


by the refraction arrival curve, e.g., the signal from ground zero. From the envelope formed 
by superimposing the arrival time curves for all signals generated by the air blast, the carliest 
arrivals are selected. This curve then defines the ground range at which outrunning occurs, 
e.g-, 920 feet in Figure 1.6. 

The measured air-blast and ground motion arrivals at 5-foot depth for Operation Tumbler 
Shot t (Reference 2) are shown in Figure 1.6. Before outrunning, the ground motion lags the 
air blast slightly since the air blast is superscismic. After outrunning, the ground motion pre- 
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cedes the calculated valucs by 2 small amount as the refracted wave travels upward to the sur- 
face. The differences are not {mportant at 5-foot depth; however, at much greater depths they 
will be significant. The 8,300-{t/sec propagation velocity 1s the result of refraction from the 
basement rock. The 1957 refraction survey did not extend sufficicntly to pick up these arrivals. 

From scismic velocitics associated with strata of known thickness, the refraction curve may 
be reconstructed. Details of this procedure may be found in any text on geological exploration. 

As the device yield increases, the near-surface seismic velocity becomes less important in 
the calculation of outrunning. The seismic velocitics at depth play the dominant role, and hence 
the overpressure at which outrunning would be expected to occur will generally increase as de- 
vice yield increases. For surface bursts over shallow weathered layers or for large yield de- 
vices detonated over deep weathered layers, the approximate values of Table 1.3 apply. This 
indication of carly outrunning for the most common materials focuses more than casual attention 
to the behavior of the ground-transmitted wave. 


1.4 SEISMIC MEASUREMENTS 


The elastic properties of a soil are related to and may be determined from an accurate knowl- 
edge of the seismic velocities. The main elastic waves which are of interest in this program 


TABLE L3 APPROXIMATE OVERPRESSURSS AT WHICH OUTRUNNING OF 
GROUND WAVE OCCURS, LARGE YIELD SURFACE BURSTS 


Formation Overpressurc Formation Overpressure 
psi psi 

Alluvium <40 Shale 650 to 2,500 

Gravel, dry 10 to 100 Limestone >1,500 

Gravel, wet 40 to 500 Metamorphic >1,000 

Sandy Clay 100 te 500 Granite >3,000 

Sandstone 500 to 2,000 


are longitudinal or compressional waves (L) and transverse or shear waves (S). Other seismic 
waves that may or may not be of interest are the Rayleigh and Love waves. 

The compressional wave (L) is characterized by an in-line particle motion with particle dis- 
placement parallel to the direction of propagation. Of the elastic waves, the compressional 
wave has the highest velocity. 

The shear wave (S) ts also linear, but the direction of particle motion-is perpendicular to the 
direction of propagation. The shear wave, because of its property of polarization, can be fur- 
ther subdivided into Sq and Sy type depending on whether the displacement is horizontal or ver- 
tical, both waves traveling with the same velocity but lower than the L wave velocity. 

The Rayleigh wave (R) is a surface wave which is a combined compressional and shear wave 
with a plane of oscillation at right angles to the surface and parallel to the direction of propaga- 
tion. Its velocity is dependent on frequency and is dispersive. Shear-wave velocity for shallow 
depths may be estimated from these surface-wave velocities. An approximate relation between 
Cg and Cp Is: 


Cr = [0.87 +0.80]Cs 


The elastic properties of a solid may be described by two elastic constants. Three of these 
constants are the bulk modulus (k), the modulus of rigidity (u), and Poisson’s ratio (9). The 
bulk modulus (Kk) is defined as the ratio between applied pressure and fractional change in volume 
fer 2 uniform hydrostatic compression. The rigidity modulus (u) is defined as the ratio between 
shear stress and shear strain. Poisson's ratio is defined as the ratio between literal contraction 
and longitudinal extension of a solid with free Lateral surfaces, It is related to {k) and (u) by 
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k+Yyu (Equation 1.1) 


The compressional and shear-wave velocities Cy, and Cs, measured by seismic methods, arc 


le 


Figure 1.7 Seismic velocity ratios versus Poisson’s ratio. 


c, -|*242]% _[_Ea-9 wa 
L P pU+0) G20) (Equation 1.2) 


_ fap =f i) w 
Cs = iB 7 [: aah] (Equation 1.3) 


Where p is the density and E Is Young’s modulus, Inserting Equations 1.2 and 1.3 in Equation 1.1 


1-2 (Co/Cp}* 


% * 3-2 (Cs/CLh 


For perfect fluids, u = 0, Cg = Qandg = °,. 
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TABLE 1.4 POISSON'S ILATIO FOR SOME MATERIALS 


A plot of ¢ versus Cy/Cy and Cy/Cy bs shown ta Figure 2.7. 


Gravel 0.47 
Locss, dry 0.44 
River Deposits 0.45 
Clay, dry 0.47 
Clay, wet 0.338 
Limestone O28 
Granite O23 
Sandstone 0.12 


For incompressible materials, k becomes large and o = %. For elastic solids, 
9<o < %. Values of o for some common materials are shown in Table 1.4. 
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Reference 7 
Reference 7 
Reference 8 


Reference 8 
Reference 7 
Reference 7 
Refcrence 7 


Chapter 2 
PROCEDURE 


Project 1.4 participated on Shot Priscilla of Operation Plumbbob. Many other projects partici- 
pated on this shot, including Project 1.3, which provided input data for this project; and Project. 
1.5, which made measurements similar to this project. There were also 2 number of structural 
projects for which Project 1.4 was designed to provide free-field Lnput data. 


2.1 PREDICTIONS 


In planning an experiment of this type, it is necessary to predict the values of the functions 
to be measured to an accuracy sufficicnt to allow the sensitivity of each channel of instrumen- 
tation to be set closely enough that satisfactory deflection may be obtained. For best resulrs 
these vajues should be within a factur of two of the true values. A greater range is acccptablic 
on channels where dual-sensitivity galvanometers are used. Predictions are also important in 
the selection of gage ratings to ensure that gages are not over-ranged, thereby introducing non- 
linearities. This section (2.1) ts presented as written before the test; it has not been modified 
in the light of the actual data obtained. 

Directly applicable data, taken In the Frenchman Flat area on underground effects of above- 
ground shots, are limited to those [rom Tumbler Shot 1 (Reference 2), and Upshot-Knothole 
Shots 9 (Reference 1) and 10 (Reference 9). The majority of these data were taken from shock 
waves at relatively low pressure levels, except some from Upshot-Knothole Shot 10 at levels 
up to 100 psi with precursor wave forms. Reference 10 gives some seismic velocities in the 
Frenchman Fiat area. . 

2.1.1 Input Predictions. Subsurface phenomena are assumed to be a result of the overpres- 
sure appearing at the surface in the immediate vicinity of the station. Table 2.1 shows some of 
the predicted parameters of interest obtained from analyses carried out on the basis of previous 
data (Reference 11). a 


2.1.2 Earth Stress. In general, vertical earth stress may be assumed to be the same as the 
applied pressure. There will be a small change in the slope of the front with depth. At Stations 
1, 2, and 3, the sharp decay of overpressure behind the front may eat away the peak slightly, 
with an estimated decrease of 25 percent, which is insignificant for range-setting purposes. 
Horizontal earth stress is assumed to one fifth the vertical. This is consistent with some of the 
previoug_experimental experience and with probable values of Poisson’s ratio. 


2.1.3 Seismic Velocity. Seismic velocity is not important per se, but it is important in pre- 
diction of strain and acceleration. If one assumes that p = 100 pcf and that o = 0.25,then with 
peak stress equal to peak surface overpressure P, Pyy (psi), Vmax = 55.6 Pyy/Cy, (in the early 
part of the wave), or Cp = 55.6 Pyy/Vinaxy where Vmax is the peak particle velocity at the 
ground surface. 

Data from Tumbler-Snapper Shot 1 and Upshot-Krothole Shot 9 in the Frenchman Flat arsa 
give average value of Cy, of 730 ft/sec at 1-foot depth, 1,220 ft/sec at 5-foot depth, and 2,460 
ft/sec at 50-foot depth. These values for Cy, at the three depths of measurement covered by the - 
datz are represented best by an cxponential form of empirical cquation: 2 


Cy = 750y°™ (Equation 2.1) 
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Where y denotes the vertical dimension. 
Reference 10 gives the following seismic velocitics in the Frenchman Flat arca: 


Depth, ft Velocity, ft/sec 
" 0 to 10 1,200 
10 to 175 2,600 
175 to 650 3,600 
Below 650 10,000 


While these values match well the average Cy results at 5S and 50 feet, they show no differ- 


ence between 1 and 5 feet. This 1s unrealistic; particularly since the abrupt velocity change at 
10 feet ts not consistent with the profile of soil observed In excavations. 


it must be observed that these data were derived from a refraction survey, where the pri- 
TABLE 2.2 SHOT PRISCILLA INPUT PREDICTIONS 


Station Ground <Arrival Horizontal Wave 


Range Time Velocity Form Pmax Pi Pz Ti Th 
ft Imsec ft/sec psi psi psi msec msec 

1 450 399 6,400 6 730 7320 — tt] —_— 
2 550 117 5,000 0 600 600 — 0 — 
3 650 —s-137 6,700 0 #45 #40 — oO — 
4 730 152 6,100 1 320 230 40 5 7.5 
Ss 850 167 5,600 1 200 1100-120 § 125 
6 1,050 205 4,300 1 100 rit] 70 3 20 
7 1,350 275 3,900 1 50 300s 40 5 20 


mary data were arrival time versus ground range from small shots. When these primary data 
are reconstructed from the reported conclusions, they [it well an empirical exponential equation, 
tz = 0.0015 x8 (Equation 2.2) 
. Where x derotes the horizontal dimension. According to Blondeau’s derivation (Reference 12), 
this would correspond to a vertical variation of velocity, 
Cy, = 1,160 y *? (Equation 2.3) 


There is a considerable difference between Equations 2.1 and 2.3, but the best compromise, 
which does serious violence to neither, appears to be: 


CL = 1,000 y °* (Equation 2.4) 


2.1.4 Earth Strain. With the same assumptions used for calculations of Cy, it may be shown 
that vertical strain, S, may be expressed as: 


S = 55.6 Py/Cy? = 5.6 x 1075 Py y~*S (Equation 2.5) 


2.1.5 Earth Acceleration, Data from Operations Tumbler-Snapper and Upshot-Knothole of 


the average values of Pay/A (ratio peak overpressure in psi to peak vertical acceleration in 
: g-units) fit well the exponential equation, 


A = 1.64 Py y7O8 (Equation 2.6) 


This is for shock-wave input at pressure ievels between 10 and 25 psi. [If the rise time were 
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constant, A might be expected to vary lnverscly as the scismic velocity, e., in the same way 
a8 Vinax: 
From Equations 2.1 and 2.5, 


Vmax = 0.056 Pay y~ °°? (Equation 2.7) 
From normal mechanies, if one assumes the rise of velocity to resemble 2 half cosinc-wave, 


A = 0.0875 Vmnax/TR (Equation 2.8) 
From Equations 2.6, 2.7 and 2.8, then 


TR = 0.0033 y (Equation 2.9) 


In other words, the steeper slope of acecleration with depth is explained by the increase of 
TR with depth. Since Tp is the duration of the rise of particle velocity, V, it should also be 
the duration of the first acceleration pulse. Data on Tp yield average values of 3.25, 7.5, and 
30 msec at 1, 5, and 50 feet, respectively, corresponding to the exponential Tp = 0.00325 y%5, 
a negligible difference from Equation 2.9. 

Equation 2.6 can be used to predict peak accelerations at Stations 1, 2, and 3. At subsequent 
Stations, however, Tp is affected by the finite rise time of the input-pressure wave. In the case 
of a precursor-type wave form, the rise time of the input wave can be quite long. If one can add 
geometrically, the ‘fluence of depth and input pressure upon rise time, (Tp = vt; + Tp); then 
one may calculate Tp for each depth and calculate A from Equations 2.4 and 2.5. 


Amax = 0.00542 Py y~*2/ TR _ _ (Equation 2.10) 
Horizontal accelerations are arbitrarily taken as one third the vertical. Although experience 
to date indicates a considerable scatter in this ratio, it appears to be always less than ane. 


2.2 INSTRUMENTATION 


2.2.1 Central Station. All central-station channels of instrumentation used on Project 1.4 
were essentially identical to those used on a number of previous projects, including Operation 
Teapot Project 1.10 (Reference 13). Wiancko 3-ke oscillators supplied carrier power to the 
transducers themselves and to modified Wiancko demodulators. The demodulated signal was 
applied to William J. Miller Corporation oscillograph recorders. Provisions were included 
for applying automatically a synthetic calibrating signal to each channel immediately prior to 
zero time to compare the final deflection on the record with the deflection produced by the same 
signal at the time of calibration. An accurate timing signal of 100 cps and 1,000 cps was also 
applied to all recorders simultaneously from a single source having a time accuracy of better 
than 10 parts per million. This provided means for time correlation of records to a high degree 
of accuracy. Since the same central station equipment was used for Projects 1.3, 1.7, and 3.5, 
this time correlation could be extended to include data obtained on Project 1.4. 

The prime power supply for all instruments during the shot was a bank of storage batteries. 
Suitable converters were used to produce 115-volts ac for components requiring this power 
source. An individual converter was used for each rectifier power supply, thus minimizing 
the possibility of gross failure due to converter failure. 

On this project, 62 gage channels were connected. Of these, 22 were connected to dual- 
recording systems, consisting of one galvanometer on each of two recorders. These dual 
channels were assigned to those gages which were considcred to be most important to minimize 
loss of important data due to any recorder failure. On 6 of these 22 channels, one of the gal- 
vanometers used had a natural frequency of 200 cps, and the remainder had a natural frequency 
of 300 cps. The channels incorporating one 200-cps galvanometer were used on gages whera 
uncertainty of the predicted peak was greatest and where the expected signal was such as not to 
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be degraded appreciably by the reduced response of the lower frequency galvanometer. Since 
there was an appreciable difference in the sensitivity of the galvanometers thus used on a single 
channe!, a wider range of input signal could be accommodated without loss of data (provided both 
recorders operated properly). 

Instruments were powercd at sultable times before zcro timc by Edgerton, Germeshausen, 
and Grier (EG&G) relay circuits, with lock-in relays controlled by a time-delay relay to continue 
operation for approximatcly 1 minute after zero time, even though EG&G relays dropped out 
sooner. Utmost attention was paid to circultry and procedures to Insure maximum reliability 
of operation. Dual-relay contacts or dual relays were used wherever feasible. 

A multipen recorder was connected to provide a record of time and scquence of operation of 
various elements so that any failure which might occur cauld be traced to its source in a posttest 
study. 

The recording shelter was buried to a depth sufficient to reduce the integrated radiation dosage 
within the shelter to below 10 r to avoid radiation fogging of the recording parer. 


2.2.2 Stress Gages. The basic earth-stress gage used in these tests was a modification of 
that originally designed by R. W. Carlson for the measurement of static stress in foundations 
and grades. The gage consisted of two flat circular plates with thin flexible edges attached at 
the edges so as to be separated by a narrow spacc fillcd with oil. A pressure gage was arranged 
to measure the pressure in this oil as a measure of the actual component of the stress in the me- 
dium in which the gage was buried. 

Two variations of this gage were used. For the low pressures, 2 Wiancko variable-reluctance 
pressure gage covered by a housing was used in an identical construction to those used on pre- 
vious projects (Reference 3). For pressures above 300 psi, a special diaphragm-type transducer, 
manufactured by Ultradyne Engineering Company, was used in place of the Wiancko transducer 
(Figure 2.1) This variable-reluctance transducer was electrically similar to the Wiancko, but 
was mechanically smaller in size so that the projecting housing was smaller, 


2.2.3 Strain Gages. The earth-strain gages used were of a new type designed for the purpose 
and shown in schematic form in Figure 2.2. The gage is shown unassembled in Figure 2.3 and is 
shown installed in Figure 2.4. Fundamentally, the gage consisted of a linear differential trans- 
former, which measured the change in spacing between two anchors, separated by 2 or 3 feet, 
set into the side of a large hole before the hole was backfilled. One of these anchors supported 
the body of the transformer enclosed in a protective housing, and the other carried a2 light tubu- 
lar rod, at the far end of which was attached the movable core of the transformer. The position 
af this core could be adjusted after installation. For calibration, the adjusting nuts could be used 
as a micrometer by turning them one full turn at a time. The linear differential transformer was 
connected through a transformer to the normal Wiancko half-bridge circuit by the standard three- 
conductor cable. 


2.2.4 accelerometers. Accelerometers used were the standard Wiancko variable-reluctance 
accelerometers used on previous tests (Reference 2). They were enclosed in protective canis- 
ters of two types: one type carried a single vertical accelerometer; the other was used where 
horizontal as well as vertical measurements were desired and carried two mutually perpendicular 
accelerometers. 


2.2.5 Instrumcnt Response. The Wiancko gage and its associated recording system was 
basically flat down to stcady-state conditions, due to its design as 2 carrier-demodulator system. 
No corrections were required therefore for its low-frequency response. The high-frequency 
response was limited either by the characteristics of the galvanometers used or by the dynamic 
characteristics of Gre transducers. The (nominal) 300-cps galvanometcrs had an undamped 
natural frequency of 315 to 340 cps and were damped to have an overshoot of approximately 74 
pereent. This corresponded to a damping factor of approximately 0.65 critical and provided a 
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Figure 2.1 Earth stress gage: 


=> 


ue 


Figure 2.2 Schematic, 


earth strain gage. 
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nominal risc time of 1.3 msec. The nominal 200-cps gulvanometers had a correspondingly 
longer rise time of approximately 1.8 msec. 

The accelerometers varied widely in sensitivity and maximum range and consequently in un- 
damped natural frequency. In general, those used at high-pressure statlons and near the surface 
were low-sensitivity, hich-range Instruments with a high natural frequency; and those used at 
low-pressure stations, particularly at greater depths, were more sensitive instruments with 
lower natural frequencies. Where the high-range instruments were used, the overall frequency 
response was limited by the galvanomcter response. Where the low-range instruments were 
used, the limiting frequency response was generally that of the accelerometer. 

The frequency-response characteristics of the Carlson stress gages were difficult to deter- 


Figure 2.3 Earth strain gage, unassembled. Figure 2.4 Earth strain gage, installed. 


mine implicitly since they were affected significantly by loading of the earth on the gages; 
however, the basic gage was known to have a response similar to that of the pressure gage alone. 
It is believed that the overall response was limited by that of the galvanometers. Similarly, the 
response time of the earth-strain gages was difficult to ascertain exactly, but measurements 
indicated that for the short span used the response time was far shorter than the rise times in- 
dicated on the final records; therefore, no distortion from this cause was attributed to the gages. 


2.2.6 Calibration. Each gage was calibrated in the field after it had been connected to its 
associated cable and recording equipment and immediately prior to its final installation in the 
earth. Carlson stress pages were calibrated by the application of direct air pressure. Accel- 
erometers were calibrated by the use of a spin table which produced accelerations up to about 
200 g. Where higher accelerations than 200 g¢ were anticipated, gages were calibrated to that 
figure in the field, and this caiibration was extrapolated on the basis of linearity checks made 
previously in the laboratory. Earth-strain gages were calibrated by the introduction of directly- 
measured deflections on the gages in place. 

In the calibration procedure several deflections ranging from zero to well above the expected 
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peak (where possible) were applied to each gage In scquence. Each galvanometer deflection was 
noted and recorded. In addition, the deficction caused by an artifictal signal (cal signal) injected 
into the gage circuit was recorded. From the former deflection, a calibration curve of deflee- 
tion versus peak reading could be constructed. The cal simnal served to correct for any changes 
of sensilivity of the recording system between calibration and the final test since an identical 
signal was injected on the iinal record about 10 seconds before zero time. 


2.3 EXPERIMENT PLAN 


2.3.1 Gage Placement. Gages were located at several depths in two sizes of holes. At five 
Stations, the holes were 30 inches in diamcter by 10 feet decp, and at two other stations the 
holes were 5 feet in diameter by 50 fect deep. Stress gages and accelerometers were placed 
in small excavations in the sides of the holes, and backfill material was carefully tamped around 
the gages by hand to ensure proper contact with the formation before the general backfilling pro- 
ceeded. Strain-gage anchors were cemented into the walls of the large holes to make good con- 
tact with the undisturbed formation and to minimize the effect of the difference between charac- 
teristics of the backfill and of the original formation. In addition, every effort was made in 
backfilling to return the material as nearly as possible ta {ts original condition. 

The specifications for preparation of backfill material and for tamping procedures were set 
up by Project 3.8 (Reference 14) on the basis of laboratory tests conducted at Waterways Exper- 
iment Station (WES) on samples of material from the area. These specifications were directed 
toward restoring the dynamic modulus rather than the density and required careful control of 
the water content at slightly below Proctor optimum, with somewhat greater than normal tamp- 
ing effort. During backfilling, control of procedure was assisted by the frequent sampling and 
analysis of the backfill in the field by Project 3.8. Samples for record were also taken by that 
project and are reported in Reference 14. 


2.3.2 Gage Coding. For identification of channels and recorded traces with their proper 
gages, a systematic coding was adopted. A station number was assigned to each gage station; 
these mambers were used as a first part of the gage code. The second part of the gage code 
was a letter indicating the type of measurement. For this project V was used for vertical 
acceleration, H for horizontal acceleration, CV for vertical stress, CH for horizontal stress, 
and SV for vertical strain. A third part of the code indicated the depth of the gage (in feet) 
below the surface. 

Typical gage code numbers would then be 3V5 [or Station 3 vertical acceleration at 5-foot 
depth; 4CV20 for Station 4 vertical stress at 30-foot depth, etc. 


2.3.3 Gage Layout. The gage layout (Figure 2.5 and Table 2.2) was selected to provide the 
Maximum of basic data on phenomenology and at the same time to provide maximum coordiration 
with other projects. Stations 4 and 6 were placed at ground ranges where predicted peaks of ap- 
plied pressure of approximately 300 and 100 psi, respectively, were expected. At these loca- 
tions, measurements were made of acceleration, stress, and strain, at a number of depths down 
to 50 feet, Including two measurements cach of horizontal accelex ation and horizontal stress. 
The remaining stations were chosen to correspond to aboveground stations of Project 1.3 and to 
cover a wide range of input pressure levels. At these stations, measurements were made of 
vertical acceleration and stress at depths of 5 and 10 fect only. 


2.3.4 Seismic Measurements. It was considered desirable to obtain data on scismic propa~ 
gatian velocities und their variation with depth, particularly in the first 100 fect, to assist in 
interpretation of final data from this project. To obtain these seismic data, a program was 
conducted in the area during the preparational phase of operations prior to the skot. Figure 2.6 
shows the location of shot holes and lines used in this work. Shot Points 1 and 2 were 200-foot 
holes, originally drilled dry, but flooded before use (Section 2.4). Shot Point J was a similar 
100-fost hole drilled iater to check the results from the first two. In these holes, small charces 
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TABLE 2.2 GAGE LAYOUT 


ppk, purls per thousand. 
Statlon Grou Gage Cude Predicted Galvanomoter Station = Qrounl Guyo Cole Predicted Galyvanometer 


ft 
4160 


v5 
§v10 
ICY5 
icVv10 


2V5 
2vlo 
2CV6 
2cvl0 
uV5 
aV10 
aCVE 
sCvlo 


dvi 
1V6 
4V10 
4V20 
4V30 
4v50 


41110 
41160 
4cvl 
4CV6 
4cvlo 
4CV20 
4cvi0 
acv50 
4CH10 
4CH50 


48V6 

43V10 
48V20 
48V30 
48V50 


270g 
140 g 
160 pul 
760 pal 


a) & 
10g 
600 pel 
600 pul 


460 pal 
460 psi 


180 g 
62g 
36g 


220 pat 
320 pal 
320 pal 
80 pal 
60 pal 


12 ppk 
9 ppk 
6.1 ppk 
6.8 ppk 

16. ppk 


Frequeno' 


400 
400 
100 
100 


400 
300 
300 
300 
300 
300 
300 
100 


300 
300 
300 
100 
300 
300 


900 


cps 


300 


300 
300 
400 


300 
300 
400 
400 
300 


$00 
300 
300 
300 
300 


200 
200 
200 


100 


300 
300 
300 


300 


300 
200 
300 


6 


6 


7 


860 


6V5 
fV10 
5CVE 
6CV10 
evi 
6V6 
WV 10 
6V20 
oVva0 
6V50 


allo 
all6o 
6cvi 
8CV6 
6cvi10 


acvz0 
6GCV30 
8UV50 
6Cillo 
8CH5O 


6SV5 

GSV10 
OSV20 
uSVI0 
g3V50 


7V6 
TV10 
IC V6 
7CV10 


39 ¥ 
223 
200 pul 
200 pel 


Gig 
22g 
iy 
Ug 
60g 
ud 9 


4.28 
Llg 
100 pul 
100 psi 
100 psi 


100 pul 
100 psl 
100 pst 
18 psl 
18 psi 


3.6 ppk 
2.8 ppk 
1.6 ppk 
1.8 ppk 
6.0 ppk 


lg 
6.3g 

60 pal 

80 pal 


Frequency 


Number Tango Pook q y Number Range Poak q y 
it eps 


300 
300 
300 
300 


300 
300 
400 
300 
300 
300 


300 
300 
300 
300 
300 


300 
400 
300 
300 
300 


300 
300 
300 
300 
300 


300 
300 
300 
300 


200 
200 
200 


100 


300 
300 
300 


300 


300 
400 
avu 


were set off at varlous depths &. the holes, and geophones were placed ncar the top of the hole 
aml at 50 and 100-fect from the hole to measure the arrival times to determine vertical compres- 
sional wave velocitics. 

Shot Point 4 was a 10-foot hole used in an experimental proyram In an effort to measure 
shear wave velocitics. Shot Points 1R, 2R, and 3R were also 10-foot holes uscd for a rcefrac- 
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rT #+350 Figure 2.5 Gage layout. 


tion profile ta determine horizontal velocities. In making the refraction profile, the geophone 
spread was left in one position and the shot point moved progressively out to 2 final offset of 
900 feet from the end of the spread. 


Results of these measurements are reported in Chapter 3. 


2.4 FIELD OPERATIONS 


Field operations on this project were concurrent with those for Projects 1.3, 3.5, and the 
instrumentation for Project 1.7 and were performed by the same personnel. A common record- 
ing shelter was used, and the data channels were intermingled. In most cases, common cable 
trenches were used. 

At the time the field crew arrived at the NTS, the only construction requirements completed 
for this project were the drilling of the 50-foot holes for gages and of two 200-foot holes for 
seismic measurements. The former were protected by local dikes from a subsequent generai 
flood, but the latter were covered with water for several days. 

The recording shelter was ready for occupancy 10 days after the crew arrived at the test 
site, and equipment was promptly installed. Cable-trenching operations were delayed so that, 
to meet schedules, it was necessary to start planting and backfilling of underground gages. 
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= al : Most of these gages were therefore calibrated with cables only partly in place and without the 
: full complement of cabics and gages conncctcd. 

Tt ts not certain whether this procederc has any effect on the accuracy of calibrations, but 
it was certainly not optimum. Calibration and planting of gages and backfilling of holes was 
completed about 8 weeks after project personnel arrived in the field. 

Seismic measurements were made 2s opportunitics arose. The two deep shot holes (200 feet) 
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Figure 2.6 Seismic locations. 


were shot after they had been flooded for several days. To check the possibility that the results 
were affected by water intrusion, another 100-[oot hole was drilled by Project 3.8 and was shot 
later. Additional seismic shooting was done on a noninterference basis. 

The records were recovered from the shelter on the afternoon of D-day. The recovery opera- 
tion required five persons to be in about a 300-mr/hr radiation field for approximately 20 minutes. 
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3.1 INSTRUMENTATION PERFORMANCE P 


Of the 64 gage channels installed on this project, 52 gave usabie records, although four 
records were incomplete due to cable breaks during the active period of recording. Of the four 
cables broken, three were associated with the gages at one station (Station 5, 850-foot ground 
range). The only apparent explanation for this behavior was that the finger cable trenches off 
the main 8-foot deep trench were tco ‘shallow at this station for adequate protection of the gage 
cables. : 

Of the 52 usable gage records, one record (6CV20) showed a peculiar wave form indicative 
of gage overload; although its value was limited, some information may be obtained from this 
record. 

Of the 12 channels producing no record, six were lost due to failure of one oscillograph to 
pall recording paper throughout the run. This camera ran for about 10 seconds after turn-on 
(until minus 5 seconds) before the recording paper tore and jammed. The other six failures 
were apparently due-to breaks in the cables or faulty cable-plug connections. Some cf these 
difficulties were detected before the shot, but too late to repair; some were apparently caused 
by strains or jarring introduced during the backfilling operations. It appears that the strain- 
gage design used on this project was particularly susceptible to electrical shorts and/or cable - 
continuity difficulties. Much of this might be eliminated by an improved cable-connection system. 

There was no definite evidence of record distarbance caused by the electromagnetic transient 
at zero time. All traces remained within the boundaries of the recording paper. Some gages, 
particularly the stress and intermediate depths, produced small deflections which reduced the 
accuracy of the data obtained. 

No evidence of radiation fogging was observed on any of the recordings obtained on this project. 

The data obtaincd for Project 1.7 was transmitted to the appropriate agency for amalysis. 


3.2 DATA REDUCTION PROCEDURES 


3.2.1 General. After each gage record was identified on the oscillograms, they were read 
(inches of deflection of record versus time) with an electromechiaunical reader, Benson-Lehner 
“Oscar” Model J. The reader output was fed into an IBM card punch, which produced the data 
cards. These deflection versus time data cards, along with appropriate calibration cards, were 
processed by an IBM Model 650 electronic computer. The final reduced data came out in the 
form of parameter (e.g., acceleration) versus time listings corresponding to each gage record. 
These listings were then plotted to give data upon which this report is based. 


3.2.2 Integration Procedurcs. Rt was desirable for the earth acceleration versus time rec- 
ords to successively integrate the results to obtain, first, the particle velocity versus time, 
and second, the particle displacement versus time. It becomes apparent after only a few at- 
tempts at this integration process that there is a good deal of judgment involved in obtaining a 
meaningful result. The main problems and their solutions are discussed fully in Appendix C of 
this report. 

‘Suffice it to say that the Integration procedure involved the following operations: 

1. Integration of the as-read acceleration-time record ta obtain as-read velocity-time. 
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2 Adjustment of the as-read velocity-time baseline to obtain zero velocity at an appropriate 
time; called corrected velocity-time. 

3. Adjustment of the accelcration-time baseline tv be consistent with the velocity-time ad- 
justment cited in [tem 2 above. : 

4. Integration of the adjusted velocity-time record to obtain the displacement-time trace, 
ealied corrected displacement-time. 

As explained in Appendix C, the corrections as applicd usually have little effect upon the as- 
read acceleration, more on the 23-read velocity, and most upon the displacement. 

The results of the Project 1.4 experiment will be presented in terms of the corrected values 
only; however, the figures and tables will indicate the magnitude of the corrections (Table 3.1). 


3.3 GAGE RECORDS AND TABLES OF RESULTS 


Figures 3.1 through 3.24 present the significant portions of the replotted gage recards ob- 
tained on this project. Figure 3.1 presents the overpressure-time records obtained on Project 
1.3 and represents the input air pressure at the ground surface on this shot. The remaining 


TABLE 3.1 SUMMARY OF PEAK SURFACE OVERPRESSURE, SHOT PRISCILLA 
CB, cable break. 


Ground Arrival Maximum Precursor Maximum Positive Fositive Wave 
Gage Range Time Pressure Time Pressure Time Phase ree Form 
Durution Impulse 
ft sec psi sec Psi sec sec psi-sec 
1B 450 0.103 None None 554.2 0.108 CB _ 0 
2B sso 0.116 39.2 0.118 366.0 0.126 CB _ 1 
3B 650 0.131 31.4 0.134 242.3 0.146 0.169 122 1 
4B 750 0.146 26.0 0.150 228.7 0.175 0.200 10.1 1 
SB 350 0.163 25.6 0.166 220.9 0.201 0.237 1L2 1 
6B 1,050 0.201 20.5 0.223 104.0 0.275 0.329 9.2 1 
7B 1,350 390.268 * 120° 0.308 * Ss9.1° 0.394° 0.783°* 6.62° 1 


* Sandia (Project L5) data. 

figures present acceleration-time, velocity-time, displacement-time, stress-time, and strain- 

time data in that order. Included in these figures are the times of arrival; designation of peak 

values; and where applicable, the air-pressure arrival (AB) at the ground surface over the gage. 
All the records are plotted to the same time scale; however, to obtain the best compromise 

between economy of space and faithful reproduction of the details of the records, it was neces- 

sary to group the records and use several different ordinate scales. 


3.3.1 QOverpressure. The pertinent overpressure-time records from Project 1.3 are pre- 
sented in Figure 3.1, and the data are listed in Table 3.1. Referring to the pressure-time plots, 
it is apparent that, although a definite precursor wave formed on the Shot Priscilla main blast- 
line, the ground level overpressure at Station 1 (450-foot range) was characterized by a clean 
(Type 0; see Appendix A) wave form. However, the longer than normal rise time of about 5 
msec gave some evidence of thermal disturbance at this station. The record at Station 2 (550- 
foot range), although incomplete due to a cable break, showed that the precursor wave formed 
betweea the 450- and 550-foot range. Subsequent pressure records documented the development 
of the precursor wave in detail. One of the most significant characteristics of this set of rec- 
ords with respect to ground motion was that the rise time to peak pressure associated with the 
main puise increased with increasing ground range and decreasing peak pressure. Also as the 
distance from ground zcro increased, the precursor wave front led the main sheck by longer 
and longer times. 


3.3.2 Earth Acceleration. Figures 3.2 through 3.6 present the replotted records of carth 
acecieration data obtained on Project 1.4. In these figures, the air-blast arrival time at cach 
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Figure 3.1--Surface overpressure versus time, Stations 1 to 7, Shot Priscilla. 
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Figure 3.2 Vertical acceleration versus time, Station 1 (GR: 450 feet); 
Station 2 (GR: 550 feet); Station 3 (GR: 650 feet); Shot Priscilla. 
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Figure 3.3 Vertical acceleration versus time, Station 4 (GR: 760 feet), Shot Priscilla. 
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Figure 3.5 Vertical acceleration versus time, Station 6 (GR: 1,050 feet), 
Shot Priscilla. 
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station is designated by the vertical line labelled AB on cach plot, and the arrival time, time of 
peak, and peak accelcration are indicated. In addition, the magnitude of the baseline correction 
applied to cach acceleration-time record is indicated by a dashed linc which designates the base- 
Une before correction. 

The vertical acceleration versus time curves are mainly characterized by 2 singJe sharp peak 
of acccicration in the downward direction, often preceded by minor disturbances. It is obvious 
that the latter are produced by the precursor and the former by the larger main peak of 2ir blast. 
The duration of the main sharp peak increased with ground range because of increased rise time 
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Figure 3.6 Horizontal acceleration versus time, Station 6 (GR: 1,050 feet); vertical 
acceleration versus time, Station 7 (GR: 1,350 feet); Shot Priscilla. 


of the input wave; it increased with depth, because of modification of the wave with travel throuzh 
the earth. At the deeper gaces, particularly at Station 6, the wave form became more complex, 
without entirely losing these characteristics. The horizontal accelerations showed somewhat 
similar wave forms, but with the first and major deflections positive (outward from ground zero). 
They tended to be slightly more oscillatory in character, with rather pronounced negative peaks 
following the major positive ones. 

Table 3.2 presents the corrected acceleration data, Ina few cases, the criteria used for 
making the integration corrections were in some doubt. Therefore, two independent choices 
were made, and both values are included in the table (see Section 3.3.3 for more details). 
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TABLE 32.2 SUMMARY OF CORRECTED ACCELERATION, VELOCITY AND DISPLACEMENT DATA, SHOT PRISCILLA 
NR, no record; (1) and (2) indicate alternate baseline choloes. 


a —tCtCsC*S™tC«C™t*C*C*C*CFGGLOG Acceleralinn ___ _Correcled Veloclty Corrected Displacement splacemeat 
Ground Gage Arrival Negativa Time of Time of Time of Time of Time of 
Gage Range Depth Time — Preunoe Precursor Ircahiie: Meximum Herm moxtmurs MOM ptaxtmum OMY testmum 
Acceleration Acceleration FO8MVE — positive NPBA Hes ative Negallve Negatlve 
t t #0 7 seo Z 00 z a0 it7aeo eee tt se¢ 
1v8 480 8 0.110 _- - 37.2 0,124 387.7 0.114 25.3 0.216 0.499 0.232 
1V1o 460 10 0,116 — - a4 0.127 =. 226.0 0,120 14.6 0.124 0.349 0.226 
2V5 650 6 0.120 _ _ 214 0.151 193.4 0,130 16,0 0.136 0.686 6.220 
2v10 660 10 0.128 _ _- 22.4 0.150 161.9 0.136 20.2 0.143 0.964 0.260 
avs 660 5 0.137 24.2 0.106 16.9 0.174 120.8 0.161 16.3 0.166 0.639 0.432 
3V10 e60 10 0.142 22.8 0.149 13.8 0.169 129.3 0.166 12.1 0.164 0.555 0.272 
4V1 760 1 0.147 40.6 0.148 26.2 0.184 168.3 0.170 16.8 0.160 0.664 0.264 
4V6 760 6 0.160 19.2 0.168 16.4 0.105 64.7 0.174 16.2 0.165 0.763 0.346 
4V10 760 10 0,166 9.39 0.162 14.0 0.198 39.2 0.180 10.9 0.183 0.621 6.324 
4V20A 780 20 0.162 _ - 14.3 0.162 46.5 0.189 7.04 0.206 0.466 6.338 
AV30A 760 40 0.170 _— - 9.70 0.220 38.1 0.197 7.02 0.216 0.363 0.324 
4V60 760 $0 0.163 _ _ 3.99 0.240 13.4 0.212 4.96 0.230 0,311 0.372 
AV50A 760 $0 0,183 - - 3.93 0.241 12.6 0.212 4.66 0.229 0.269 0.372 
410 7180 rt] 0,152 _— _ 14.3 0.178 6,78 6.194 1.42 0.184 0.076 0.876 
41150 760 60 0.186 _ _ 6.07 0.212 2.68 0.222 0,985 0.216 0.014 0.232 
6V5* 880 6 0.170 16.2 Q174 8.37 0.218 24.3 0,208 10.6 0.214 _ _— 
6V10¢ 850 10 0.176 6.88 0.184 4,66 0.196 16.7 0,208 6.47 0.220 - — 
GVEQi) 1,050 1 0.204 16.0 w 0,206 7.61 0,284 21.8 0.272 6.12 0.274 0.218 0.343 
{2} 16.2 0,206 9.42 0,284 21.7 0.272 6.68 0.274 0,269 0.404 
6V5(1) 1,080 6 0.207 10.4 0.2142 7.01 0,290 16.7 0.274 6.02 0,262 0.287 0.986 
(2) 16,6 0.212 6.64 0,200 16.8 0.274 6.18 0.262 0.307 . 0.418 
6V10 1,080 10 NR _ _ - - _ _ - — — - 
6V20(1) 1,050 20° 0,217 4,86 0.224 2.62 0.315 10.5 0.277 4,07 0.268 0.308 0.814 
(2) 4.87 0,224 2,52 0.316 10.8 * 0,277 4.07 0.299 0.307 0.436 
6V20A 1,080 20 0.217 4,46 0,226 2.66 0,318 9.63 0,278 3.98 0.208 0.276 0.586 
(2) 4.82 0,226 2.62 0.316 . 10.0 0.278 4.08 0.300 0.293 0.398 
6VI30 1,060 10 0.220 1,64 0.228 6.08 0.314 6.63 0.286 3.68 0.312 0.214 0.672 
6V30A 1,080 30 0.220 1.64 0.220 7.84 0.314 8.40 0,288 3,66 0,312 0.209 0.622 
6V60(1) 1,080 60 0.292 1,69 0.242 2.71 0.337 2.72 0,296 2.30 0.932 0,169 0.426 
{2} 1.60 0.242 2.69 0.337 2.74 0.208 2,36 0,332 0.182 0.462 
6H10 1,050 10 0.212 _ - 1.31 0.274 2.80 0.299 0,400 0.266 0,008 0.324 
6HSO 1,050 60 0,210 _— _ 2.84 0.296 1.65 0.322 0,920 9.303 6.022 0.826 
qV5 1,950 6 0,269 4.03 0.278 4.40 0.403 9.16 0.304 2.84 0.400 0,239 0.646 
TV10(1} 1,350 10 0,276 2.49 0,263 1.49 0.231 4.64 0.406 1.49 0.405 0,090 0.628 
{2) 2,080 2.44 0.283 1.68 0.201 4.78 0.400 1,67 0.406 0.164 0.662 


* Uncorrected values (cable break at 0.2200 ao). 


Where dual channels were employed on a single gage, the two results are shown in the table; 
however, only onc of the pair ts plotted in the figures. ' 


3.3.3 Earth Velocity. The corrected earth velocity versus time plots are presented in Fig- 
ures 3.7 through 3.13. The magnitude of the bascline change necessary to obtain zero velocity 
at a specificd time is indicated on each curve; the dashed linc designates the baseline of the as- 
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Figure 3.7 Vertical velocity versus time, Siation 1(GR: 450 feet), Shot Priscilla. 


read first integration. In general, the wave forms of velocity-time are similar. At the close- 
in stations, where the precursor is just forming, the curves rapidly rise to maximum velocity. 
However, as the precursor develeps, its influence is evident on the velocity response. Also, 
the velocity peaks tend to become less sharp with increasing depth of measurement. 

The criterion for specification of the time that the velocity is zero was difficult to establish. 
After considerable study, It was decided tmt, for local air-blast induced effects only, the end 
of the overpressure positive phase would constitute a reasonable criterion for velocity equal to 
zero. Of course, this criterion necessarily was only indicative of the corrections to be made; 
each choice of time of zero velocity had te be made, taking into account any peculiarities which 
appeared on individual records. 

Ina few cases (for example: 6V5, Figure 3.12), two different bascline corrections were 
thought to be equally valid. For these, both valucs are carried through, and the two displace- 
ment values are determined. 

The data cbtained from the corrected carth velocity plots are tabulated in Table 3.2. 
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Figure 3.9 Vertical velocity versus time, Station 3(GR: 660 feet), Shot Priscilla. 


Figure 3.8 Vertical velocity versus (ime, Station 2 (GR: 560 fect), Shot Priscilla, 
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Figure 3.10 Vertical velocity versus time, Station 4 (GR: 750 feet), 
Shot Priscilla. 
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Figure 3.11 Horizontal velocity versus time, Station 4 (GR: 750 feet); 
vertical velocity versus time, Station 5 (GR: 850 feet); Shot Priscilla. 


45 ° 


SECRET 


Preseore, gal 


vey ne Sacre 


Voloalip, Ci fens 


ry me 


Figure 3,12 Vertical velocity versus time, Station 6 (GR: 1,050 feet), 
Shot Priscilla. 
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“. 3.3.4 Earth Displacement. Figures 3.14 through 3.17 present the corrected earth displice- 
Ment versus time plots, which arc obtained from the double Integration of the accelerometer 
records; the pertinent data are tabulated tn Table 3.2. Included in Table C.1 (Appendix C) is 
the approximate uncorrectcd peak displacement; that Is, the displacement obtained by double 
intexration of as-read acceleration-time without correction to zcro velocity. 

tn general, the wave forms of the displacement plots are similar, exhibiting, Ike the veloc- 
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Figure 3.13 Horizontal velocity versus time, Station 6 (GR: 1,050 fect); vertical 
velocity versus time, Station 7 (GR: 1,350 feet); Shot Priscilla. 


ity results, increasing time of rise to peak amplitude for increasing ground range. The effect 
of the precursor upon displacement is obvious at Station 6 but to a much lesser degree at the 
closer gave stations. 
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Figure 3.14 Vertical displacement versus time, Station 1 (GR: 450 feet); 
Station 2 (GR: 550 feet); Station 3 (GR: 650 feet); Shot Priscilla. 
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Figure 3.15 Vertical and horizontal displacement versus time, Station 4 
(GR: 750 feet); Station 5 (GR: 850 feet); Shot Priscilla. 
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Figure 3.16 Vertical displacement versus time, Station 6 (GR: 1,050 feet), 
Shot Priscilla. 
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Figure 3.17 Horizontal displacement versus time, station 6 (GR: 1,050 feet); 
vertical displacement versus time, Station 7 (GR: 1,350 feet); Shot Priscilla. 
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3.3.5 Earth Stress. The carth stress-time plots from Project 1.4 are shown in Figures 
3.18 through 3.23. Included in cach figure, for comparison, is a plot of the overpressure-time 
record obtained at each station. 

Tt {3 seen that, with few exceptions, the wave forms of these plots, particularly those from 
Shallow gages, are similar to those of applied air-pressure. At the grcater depths, the time of 
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Figure 3.18 Vertical stress versus time, Station 1(GR: 450 feet); 
Station 2 (GR: 550 feet); Shot Priscilla. 


rise to peax stress is increased, and at Station 6 the effect of the precursor wave is evidetu. 
Three records, 6CVSB, 6GCV!0B, and 6CV20B, obtained from Project 1.7 and buried near Proj- 
ect L4, Station 6, (Figure 2.5), are included in Figure 3.22. These records will be considered 
auxiliary data and will be included in the analysis of results. It is apparent that, whereas at 
the 5-foot depth the primary and auxiliary stress records are comparable in wave form, the 
records at 10- and 20-foot depths are not similar inform. Actually, the auxiliary records 
seem more reasonable; this will be discusscd more fully in Chapter 4. 

The horizontal stress-time plots appear to exhibit the same general wave form as the vertical 
component. 

Stress data are given in Table 3.3. 


3.3.6 Earth Strain. Figure 3.24 presents the plots of the reduced earth strain records ob- 
tained, with peak values marked similarly to those of acceleration. It is seen that, while only 
two of these records are complete, three others give fragmentary information. The intermittent 
nature of these records is apparently caused by momentary contact of broken cables or connec- 
tors caused by earth motions. 

The two complete strain records show wave forms which rise to 2 peak about the time of the 
peak stress, but decay less rapidly and to a constant value which represents the permanent 
(residual) strain. This wave form is appzrently followed by the [ragmentary record of 6SV30, 
but that of 4SV30 implies a final value of strain which is negative, unless that portion of the 
record is displaced. The abrupt chanyes in strain observed on the 4SV10 plot (near 0.180 sec- 
ond) and GSV30A (near 0.280 second) make these results somewmt suspect. One would not ex- 
Peet the Strain In soil to Vary in that munner. 
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Figure 3.19 Vertical stress veraus time, 
Station 3 (GR: 650 feet), Shot Prisallla, 


Figure 3.20 Vertical strees yersus time, 
Station 4 (GR: 760 feet), Shot Prisollla, 
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TABLE 5.3 SUMMARY OF STRESS DATA, SIIOT PRISCILLA 


NR, no record; DI, duta indeterminate, und CB, cyble break. 


Gage Gfound Gage Arrival lina Time of First Maximum ices Maximum Ecler 
Range Depth Time Peak Positive Peak Positive Positive Negative Negative 
ft ft sec psi sec psi sec psi sec 
1¢V5S 450 $ 0.124 «6104.1 0.129 213.7 0.147 _ _ 
1CV10 450 10 NR NR NR NR NR NR NR 
2cvs §50 5 NR NR NR NR NR NA NR 
2Cv10 550 10 0.130 203.5 0.136 203.5 0.136 208.2 0.143 
3cvs 650 5 0.133 24.2 0.141 220.8 0.158 —_ _ 
3cvio 650 10 0.149 205.8 0.168 211.8 6.174 7.79 0.285 
4CV1 750 1 0.146 18.9 0.150 214.9 0.175 19.3 0.415 
4cvs 750 5 0.151 42.3 Q.159 181.5 0.183 34.1 0.300 
4acCvio 750 10 0.164 11.8 0.172 46.6 0.197 _— _ 
4CV10A 750 10 0.164 7.11 0.167 42.2 0.150 _ _ 
4cv20 750 20 NR NR NR NR NR NR NR 
4CV20A 750 20 0.165 28.4 0.180 28.4 0.180 _ _ 
4Cv30 750 30 NR NR NR NR NR NR NR 
4CV30A 7S0 30 NR NR NR NR NR NR NR 
4cvso 750 50 0.197 114.4 0.211 132.2 0.250 —_ _ 
4CVS50A 750 so xR NR NR NR NR NR NR 
4CH10 750 lo 0.157 20-5 0.165 94.2 0.192 13.0 0.375 
4CH50 750 50 9.199 10.0 0.212 14.9 0.247 3.98 0.207 
4CHS0A 750 50 0.199 10.4 0.220 14.7 0.248 4.36 0.207 
scv5 850 5 NR NR NR NR NR NR NR 
scvl1sc 850 le 0.181 48.6 CB 0.219 CB 48.6 CB 0.219 CB cB cB 
6cv1 1,050 1 0.204 27.6 0.227 74.5 0.273 8.97 6.515 
6cvs 1,050 3 0.207 14.8 0.228 94.4 0.281 8.58 0.530 
6CVSB 1,050 $s 0.205 133 C.222 56.5 0.279 12.2 0.445 
6écvl10 1,050 10 0.213 «DI DI $4.5 0.283 33.8 0.220 
6CV1OA 1,050 10 NR NR NR NR NR NR NR 
6CV10B 1,050 10 0.208 30.1 0.238 160.1 0.288 _— _ 
6cv20 1,059 20 NR NR NR NR NR NR NR 
GCV20A = 1,050 20 0.219 DI DI DI DI 9L3 0.260 
6CV20B 1,050 20 0.216 5.38 0.254 19.4 0.323 1.24 0.525 
6CV30 1,050 30 0.230 3.49 0.271 TAT 0.303 4.79 0.256 
6CV30A = 1,050 30 0.230 2.64 0.272 5.93 0.303 5.99 0.255 
6CVS0 1,050 Ft] 0.256 2.99 0.270 21.7 0.314 2.25 0.500 
6CVSOA 1,050 50 0.256 DI Do 23.6 0.316 2.72 0.525 
6CHL1O 1,050 10 0.216 6.44 CB 0.286 CB 6.44CB 0.286C3 cs cB 
6CH50 1,050 S50 0.281 9.640 0.298 9.64 0.298 0.438 0.485 
6CHSOA ——iL,0S0 50 0.281 8.68 0.237 8.68 0.297 _ _ 
7CVS 1,350 Ss 0.270 A.51 0.334 28.3 0.406 0.501 0.362 
7CV1LO 1,350 10 NR NR NR SR NR NR NR 
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Table 3.4 presents the basic strain data obtained. For convenience, peak alr pressure and 
peak measured stress are included. 
3.4 SEISMIC MEASUREMENTS 


The seismic measurcments were made to obtain data to aid in the interpretation of the below- 
ground effects on this project. The first phase of the program was the determination of the 
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Figure 3.21 Horizontal stress versus time, Station 4 (GR: 750 feet); vertical 
stress versus time, Station 5 (GR: 850 feet); Shot Priscilla. 


variation of velocity (Cz,) with depth, particularly in the first 100 feet. Vertically drilled holes 
were used for this purpose. Small dynamite charges were fired at various depths in these holes, 
and geophones were pliced near the top of the hole and at distances of 50 and 100 feet from the 
hole. The arrival times were measured and plotted against depth to determine the compressional 
wave velocities (Cy) in a vertical direction. 

The secoad phase of the scismic program consisted of muking a Seismic refraction profile to 
determine the seismic horizontal velocities and their variation with depth. Figure 2.6 shows the 
location of shot holes and refraction spread used in the seismic program and their relation to the 
blast line. During this part of the program an cffort was also madc to gencrate shear waves and 
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Figure 3.22 Vertical stress versus time, Station 6 (GR: 1,050 feet), Shot Priscilla. 
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Figure 3.22 Continued. 
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Figure 3.23 Vertical and horizontal stress versus time, Station 6 (GR: 1,050 feet), 7 
Shot Priscilla. 
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Figure 3.24 Vertical strain versus time, Station 4 (GR: 750 feet); 
Station 6 (GR: 1,050 feet); Shot Priscilla. 
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Figure 3.25 Travel-time data, Shot Paint 1, Frenchman Fiat, 
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Figure 3.26 Travel-time data, Shot Point 2, Frenchman Flat. 


measure shear-wave velocities; however, no success was met in the attempt to measure shear- 
wave velocity directly. However, some short spreads, using only a blasting cap and hammer 
blows for energy, werc set out and Rayleigh waves were recorded. From these Rayletgh veloc- 
ites an approximation of the surface shear-wave velocity was made. 


3.4.1 In Situ Seismic Velocitics. The seismic velocities in a vertical direction were meas- 
ured by detonating smal! dynamite charges ('/ to 154 pounds) at 2 number of depths in three shot 
holes (see Section 2.3.4) with detectors placed at the surface. The horizontal velocities were 


TABLE 3.4 SUMMARY OF STRAIN DATA, SHOT PRISCILLA 
NR, no record; and CB, cable break. 
Ground Gage Arrival Maximum 7!  ypanum Timeof  pecidar Time of 


Gage Maximum Maximum Residual 
Range Depth Time Positive cositive NeSUIve Negative FO strain 

fe ft sec Pppkr sec ppk* sec ppk* sec 

45v5 750 3 0.150 $82 9.194 _ _ 6.60 0.309 CB 

4SVSA 750 $s 0.150 8.73 0.195 — _ 7.54 0.309 CB 

4SV10 750 to 0.172 18.1 0.187 0.233 0.176 15.1 0.287 CB 
4SV20 750 20 NR NR NR NR NR NR NR 
SSV20A 750 20 NR NR NR NR NR NR NR 
4SV30 7) 30 NR NR SR SR NR NR NR 
4SV50 750 sO 0.185 7.58 0.248 —_ —_ 3.65 7.50 
4SVSOA 750 So 0.185 7.386 0.249 —_ _ 3.70 4.80 
6svs5 1,080 5 NR NR NR NR NR NR NR 
ésvio 1,050 10 NR NR SR NR NR NR NR 
€sv20 1,050 20 NR NR NR NR NR NR NR 
6SV20A 1,050 20 NR NR NR XR NR NR NR 
6SVi0 1,050 30 0.220 $.73 0.304 _— _ 0.758 7.50 
6SVIOA 1,050 30 0.220 5.99 0.304 _— —_ 0.730 4.80 
6Svso 1,050 S50 0.209 2.99 0.330 0.034 0.213 0.457 S.70 
6SVS0A 1,050 50 NR NR NR NR NR NR NR 


* Parts per thousand. 


determined from the refraction profile by setting out a line of geophones and shooting small 
charges in shallow holes at various distances in line with the geophone spread. 

The results of the vertical seismic profiles shot at the three different positions show the 
variations of the velocities with depth. Figures 3.25 through 3.27 are plots of the time-depth 
curves at the three hole positions. While there is some 20 percent lateral variation, all the 
time-depth curves are similar. (No effects of water intrusion on the vertical velocities was 
apparent in the third hole.) The average-velocity depth and the interval velocities, determined 
by the increments of the time-depth curves, are also shown in these figures. It is not implied 
that the interval velocities are as abrupt as shown, but this presentation is conventional, and 
attempts to smocth the curve are not justified by data of this type. 

The results of the seismic survey do not show the gradual increase of velocity with depth 
that was expected. The low-velocity layer in the vicinity of 40- to 50-foot depth was not indi- 
cated from previous data takeza by earlier refraction surveys (Reference 10). The refraction- 
survey method has the inherent limitation that it cannot detect a low-velocity layer which under- 
lies a high-velocity layer. When the 50-foot instrument hole was examined, there was no visual 
indication that the playa in the 40- to 50-foot low-velocity region was any different than the upper 
high-velocity material. One explanation for this low-velocity material would &c that rate of dep- 
osition for this region was faster than for the rest of the overlaying material. This higher rate 
would produce material that was not as compact and of a lower modulus and would thus have a 
lower velocity. 

The unexpected results of the vertical secismic-velocity surveys were sufficiently different 
from those assumcd in the original predictions to cause Some concern as to the validity of the 
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predictions. There was insufficient backyround ay to the cfiects that the velocity inversion 

would have on the stress and acceleration to justify thelr modification. The strain, however, 
was considered to be more directly affected by the local modulus, and these predictions were 
modified to take the selsmic results into account. The revised values are shown in Table 2.1. 

The results of the refraction profile show four layers having different horizontal velocities 
to a depth of 200 feet. The surface layer to a depth of 20 feet showed an average velacity of 
1,050 + 100 ft/sec, the second layer, which was 42 feet thick, had an average velocity of 2,260 
ft/sec and the third layer, 106 fect thick, had an average velocity of 2,665 {t/sec. The fourth 
layer was penetrated only 2 few fcet and showed a velocity of 3,340 ft/sec. If the scismlc pro- 
file had been extended beyond 900 feet along the surface, the depth and velocity of deeper layers 
could have been determined. 

Figure 3.28 shows the travel curve of the refraction profile and 2 wave-front diagram indi- 
cating the depths and velocities of the various refraction horizons. The aorizontal velocities 
shown by the scismic-rcfraction work < re somewhat greater than the vertical velocities. This 
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Figure 3.2T Travel-time data, Shot Point 3, Frenchman Flat. 
is normal and to be expected in multilayered formations. As stated before, a break in the time- 
travel curve occurs only when a higher velocity layer is below a layer of low velocity. The ver- 
tical velocities probably give better overall velocity data than the refraction velocities which 
only indicate velocities in specific layers. 

A close examination was made of the refraction records to see if any shear waves were re- 
corded. No definite arrivals were observed that could be identified as shear waves. However, 
two short spreads were recorded using geophone intervals of only 20 feet and hammer blows 
and a blasting cap as a source of energy, and Rayleigh waves were observed on these short re- 
fraction profiles. The results of this survey are shown in Figure 3.29, which shows a Rayleigh 
wave velocity of 575 ft/sec + 15 percent and a Cy, of 1,520 ft/sec + 15 percent. This gave a 
ratio of Cp/Cy, = 0.378 (Figure 3.29) which gives a Poisson’s ratio (a) of 0.40 for the surface 
material to a depth of approximately 15 feet; esing the extreme values of velocity, one obtains 
0.22 = o = 0.44. 


3.4.2 Velocities in Backfill. During the backfilling of the deep guge holes at Stations 4 and 
6, a geophone was planted in each hole ncar the wall at a depth of 50 fect. Late in the preshot 
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operation after backfilling was compicted, measurements were made of the travel time of the 

seismic Lmpulse from the geophone to the surface. ‘ 
Small charges consisting of only 2 blasting cap were fired on the surface. One charge was : 

shot on the opposite side of the hole from that where the geophone was plantcd and the other 

charge 5 feet away from the edge of the hole above the gcophone (also 5 fcet off vertical). The 
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Figure 3.29 Refraction survey, short-span instrument array, Frenchman Fiat. 


travel times through the undisturbed material were 33 and 36 msec and through the backfilled 
material they were 36 and 37 msec. These travel times correspond to average velocities of 
1,390 + 40 ft/sec and 1,450 + 60 ft/sec. The average travel time for the first 50 feet from Shot 
Points 1, 2, and 3 showed a time of 24.2 msec which corresponded to an average velocity of 
2,060 ft/sec. These data imply that the backfill failed to reach the original modulus and that 
the effect of backfilling and moisture content also reached somewhat beyond the edge of the hole. 
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Chapter 4 
DISCUSSION 


After a brief description of acceleration and velocity wave form types, the Shot Priscilla ground 
shock data will be discussed. Where applicable, the results of the Tumbler, Upshot-Knathole, 
and Jangle shots (Appendix B) will be compared with those from Shot Priscilla. The main topics 
to be discussed are attenuation of ground shock with depth, carth response as affected by the 
character of air-blast pressure input, local effects versus effects from remote sources, and 


response spectra. Also included will be a brief analysis of deduced stress-strain relations for 
soll. 


4.1 ACCELERATION AND VELOCITY WAVE FORMS 


The most typical characteristic of ground accelerations induced by the air blast is that there 
ts no single ideal wave form. In the elastic model, the accelerations caused by the direct wave 
appear proportional to the time rate of change of overpressure; hence, one would expect the 
character of the accelerations to change with overpressure wave form. The integral of vertical 
acceleration, vertical particle velocity, is perhzps the most familiar of all ground-motion pa- 
rameters Since it bears a direct relationship to the overpressure. Six of the more predominant 
acceleration and velocity wave form types characteristic of the NTS are shown schematically in 
Figure 4.1. Examples of these wave form types taken from various tests are listed in Table 4.1. 
When the wave. form is ideat and traveling superseismically, the downward acceleration is large 
compared with the following upward accelerations (Figure 4.1a), examples: Records 1V5, 1V10, 
and 2V5 (Figure 3.2). Tke resulting velocity is similar to the air pressure, falling off somewhat 
more rapidly than the pressure and becoming zero before the end of the positive phase of the air 
pressure. At lower peak overpressures, the first downward acceleration is followed by an up- 
ward peak acceleration of nearly equal or sometimes greater magnitude (Figure 4.1b). This 
type of wave form is common in small high-explosives charge work even at high overpressures. 

The resulting velocity decreases more rapidly than in Type a and becomes positive prior to 
the end of the positive phase. The peak upward velocity is always much smaller than the peak 
downward velocity, however. The first indication that the air blast is traveling transseismically 
is an outrunning of the downward acceleration (Figure 4.1c) compared with the sharp onset of 
acceleration illustrated previously, example: Record 4V20A (Figure 3.3). This type of outrun- 
ning, although infrequently observed during nuclear tests, is common on small high-explosives 
shots. 

The character of the greund acceleration in the precursor region changes markedly as the 
precursor develops and decays. Figure 4.1d illustrates a wave form typical of the early stages 
of precursor development, examples: 4V1, 4V5, and 5V5 (Figures 3.3 and 3.4). The precursor 
induces accelerations similar to those of Figure 4.1b. The main pressure wave induces larger 
(since the pressure is much greater) downward accelerations. The peak upward accelerations 
are low since the pressures are large in this region. As the precursor transgresses its various 
Stages of development, the air pressure wave form becomes violent, producing the high frequen- 
cy acceleration record of Figure 4.le, examples: 6V1, 6V5 (Figure 3.5). As the blast wave 
cleans up, these accelerations may become small in magnitude. 

When the refracted ground wave outruns the local wave, the first peak vertical acceleration 
ts positive or upward as in Figure 4.1f. The frequencies associated with the remote source 
accelerations are much lower than those of the local effect; onset of acecleration is gradual and 
in many cases hardly perceptible; and successive peaks grow in magnitude and then decay. The 
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signals often last for lwo or three times the positive phase duration of the air blast wave. Al- 
though the accelerations from remote sources are gencrally small compared with those of the 
local wave, their low frequency results in velocities and displacements which are comparabic 
with those of the local wave. 

The arrival of the local effect is tdcntified by a sudden increase in velocity (Figure 4.2), the 
velocity jump, as it is called in British weapons effects reports. Since the velocity jump may 
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Figure 4.1 Schematic diagrams of vertical acceleration and velocity wave forms. 


had 


be superimposed on ascending or descending portions of the remote effect, it represents the 
most probable peak downward value of the velocity. Furthermore, it bears a direct relationship 
to the overpressure in the same manner as the peak downward velocity tn the superseismic case. 
The influence of the refracted wave not only appears 2head of but alse behind the arrival of the 
local wave. Because of this, the peak upward velocities become comparable with the velocity 
jump, Figure 4.3. 

Even though the blast wave may be superseismic, f.e., outrunning is not observed, the sig- 
nal produced oy the refracted wave may be sufficient to modify the velocities in the later per- 
tions of the acceleration pulse. Thus, the remote source wave can have marked effects on the 
peak displacements, effects which are not observed in the peak velocity. 

At depth, outrunning occurs carlier than at the surface, and the accelerations of the direct 
wave 2re attenuated. In such cases, the acceleration wave type becomes extremely confused, 
particularly in the precursor rcgion where types d, f, and possibly c emerge. This combina- 
tion may occur at depths as small as 20 fect, cf., Shot Priscilla gages 4V20A and 6V30A (Fig- 
ures 3.3 and 3.5). 


4.2 SIGNALS FROM REMOTE SOURCES 


The origin of sigrals from sources other than the local air-blast slap has becn discussed 
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for the general case in Section 1.3. Figure 4.4 shows the air pressure time of arrival plot for 
the Priscilla event; also pluited on the figure ure the data from the refractlon survey taken at 
Frenchman Flat as a part of Project 1.4. Sliding the setsmic curve up the air-blast curve, one 
finds that the surface level outrunning ground range for Priscilla was about 2,500 feet, due to a 
signal originating approximately 600 feet closer to ground zero. However, the plot of Figure 
4.4 does not tell the whole story; !t is possible for: (1) outrunning to occur at closer ranges for 
mcasurements taken underground; and (2) refracted signals to be manifest on records after the 
arrival of the first signal. 

Therefore, although the Priscilla acceleration results (measurements out to 1,350-foot ground 
range) show no outrunning in the strict sense of the term, they show that refracted signals from 
remote sources were present, particularly {rom the deeper gages. Notable examples are found 


TABLE 4.1 EXAMPLES OF VERTICAL ACCELERATION WAVE 
FORM TYPES SHOWN IN FIGURE 4.1 


Wave Form Example * Shot Overpressure 
psi 

2 Ivs Plumbbob, Priscilla 5534 
14Vi Upshot-Knothole 9 21.5 

b 1TV1 Upshot-Knothole 9 1L5 
5v Mole Round 104 78 

ce 4v Mole Round 208 >160 

d 4V1 Plumbbob, Priscilla 104 
qvs Plumbbob, Priscilla 33 
svi Upshot-Knothole 10 300 
Iv Tumbler 4 4s 

e 1vTv1 Upshot-Knothole 10 14.5 
ovi Upshot-Knaothbole 10 8.1 

£ 8v Tumbler 1 §.2 
sv Tumbler 2 3.4 
9v2 Mole Round 308 13.8 


* Gage designation corresponds to that given in reports of these shots. 


at Stations 6 and 7 (Figures 3.5 and 3.6). Record 6V30A exhibits the type of behavior expected 
when the signals from remote sources arrive after the local precursor slap and the main wave 
slap; however, when this occurs, it is almost impossible to identify the effect of the refracted 
wave. Similar behavior is evident on the TV5 and 7V10 records (Figure 3.6); it is likely that the 
low frequency, small amplitude acceleration oscillations beyond about 0.5 second on these rec- 
ords are due to remote source signals. 

tt is readily apparent that the remote source disturbances are extremely bothersome when 
one is integrating acceleration-time records to obtain velocity and displacement versus time. 
Any criteria which are ordinarily useful for determining the time of zero velocity (Appendix C), 
may be rendered void by the presence of signals from 2 remote source. 


4.3 ATTENUATION OF GROUND SHOCK WITH DEPTH 


By far the major portion of data which may be used to determine the attenuation of ground 
motion with depth are the measurements on Shot Priscilla. For this reason, observations made 
on other shots are discussed concurrently with the Shot Priscilla data. 


4.3.1 Acceleration (Attenuation with Depth). The attenuation with depth of maximum down- 
ward acceleration measured on Priscilla is summarized in Figure 4.5, Both the Stanford Re- 
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Figure 4.2 Uncorrected and corrected vertical velocities, 5-psi overpressure, 
Tumbler Shot 1 (arrow shows air-blast arrival). 
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Figure 4.3 Vertical velocity versus ground range, Tumbler Shot 1, 5-foot depth 
(numbers adjacent to data are overpressures in psi). 
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search Institute (Project 1.4) and the Sandia (Project 1.5: Reference 15) data are plotted on the 
figure, Reference to Figure 4.5 shows that between 5- and 10-foot depths the attenuation varics 
between about 30 to 45 percent, except at the 550- and 650-foot ground ranges where the change 
in peak acceleration with depth Is negligible. There is no apparent explanation [or this behavior; 
it is significant that the Sandia data Indicated rather sharp attenuaticn at greater depths at a 
range of 650 feet. 

Taking the Stanford Research Institute (SRI and Sandia data together at 750- and 850-foot 
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Figure 4.4 Air-blast time of arrival, Shot Priscilla, and seismic soil survey, 
Frenchman Flat. 


ranges, it is evident that the tendency ts toward a decided decrease in acceleration attenuation 
between 10- and 30-foot depths. Wave theory gives the following formula for the particle veloc- 
ity at the interface between two materials after passage of a step pulse: 


¥ = oe ae 
i“ 
i+ oC; (Equation 4.1) 


Where @ is the stress change across the approaching wave front; p,C, is the impedarice of the 
material in front of the interface; and p,C, is the impedance beyond the interface. The stcp 
velocity across the approaching front is V = @/p,C,. Therefore, the theory suggests the fol- 
lowing: (1) If an aceclerometer is located at or near the transition from a hard to 2 soft mate- 
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rial, the peak acceleration ts increased. (2) If an accelcrometer is located at or near the 
transition from a soft to a hard material, the peak acceleration is dcercascd. 

Reference to the selsmic profile measured at Frenchman Fiat prior to Priscilla (Scction 
3.4) indicates that the prescnec of an underlying soft (lower seismic velocity) layer may explain 
the behavior at the 750- and 850-fvot ranges on Priscilla. 

Because the 10-fcot deep gage at 1,050-foat ground range was lost, there are no data between 
5 and 20 fect at this station. At 1,350-foot range, the SRI and Sandia data agree well and Indi- 
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Figure 4.5 Earth acceleration versus depth, Shot Priscilla. 


cate no pronounced change in logarithmic attenuation with depth for peak downward acceleration. 

To summarize the attenuation of peak outward horizontal acceleration between 10- and 50- 
foot depths, it is apparent that at the two stations instrumented (759 and 1,050 feet) the attenua- 
tion is less than for the corresponding peak downward component; numerically, it is about 40 
percent for the outward accvicration and about 80 percent for the downward. 

Peak negative accelerations on Prisciila and Jangle S, normalized against the peak accelera- 
tion at 10-foot depth, are shown in Sigure 4.6. The.10-foot depth was chosen as a basis of com- 
purison since the 3-foot depth was not common to these obscrvations. The dashed lines shown 
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in the figure are estimates of the confidence limits of the correlation, The Identical data are 
also shown In Flgure 4.7 on an exponential plot. The power law plot of Figure 4.6 appcars most 
satisfactory. 

The accelerations at 5-foot depth are 1.3 + 0.2 times those at 10-foot depth. Using this figure, 
the Jangle S and Sandia Priscilla measurements were adjusted to conform to a normalization 
against the acceleration at 5-foot depth, Figure 4.8. To complete the acceleration data summary, 
the Tumbler observations have also been plotted. The same confidence limits are shown as in 
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Figure 4.6 Attenuation of vertical acceleration with depth, 
Nevada Test Site (power law attenuation). 


Figure 4.6 and are seen to include 93 percent of al! acccleration measurements at the NTS. 

Because of the wide variation in accelcration-time wave forms, it is well known that analysis 
of peak accelerations alone has only limited value. Changes of response with depth will now be 
discussed using the velocity and displacement data. 


4,3.2 Velocity (Attenuation with Depth). In Figure 4.9 are shown the vertical velocity jump 
at depth referred to the velocity jump at 5-foot depth. This manner of presentation has been 
chosen to tie into the correlation of Figure 4.27. An exponential decay lw appears to represent 
the data better than a power law decay. There is nothing fundamental 2bnut this choice of coor- 
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dinates; the choice was simply a matter of convenicnce in providing a straight linc <xtrapolation 
to greater depths. 

Also shown in Figure 4.9 are five datum points from Tumbler. Threc of these observations 
are in excellent agrcement with Priscilla data; two are 50 percent larger. 

On Priscilla, at the two stations (peak overpressures of 275 and 100 psi) where measure- 
ments were obtained down to 50-foot depth, attenuation of penk downward velocity is similar 
for these overpressure inputs. There is some indication on the figure that the underlying soft 
layer affected the peak downward velocity at the 100-psi station; however, at 275 psi, the effect 
observed on the peak acceleration (Figure 4.5) is almost completely eliminated. 

The peak outward velocity data from Priscilla shows somewhat less attenuation with depth 
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Figure 4.7 Attenuation of vertical acceleration with depth, 
Nevada Test Site (exponential law attenuation). 


than the downward component at 275-psi overpressure; at the larger ground range (100-psi 
overpressure), the outward velocity at 50-foot depth ts actually more than double that measured 
at 10 feet. This Latter effcet can be explained on the basis of the intlucnce of ground motion sig- 
nals from remote sources closer to ground zero. 


The vertical velocity-time curve is similar to the overpressure-time curve when the air- 
blast wave is superscismic. The velocity peak is, however, more rounded and progressively 
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Figure 4.8 Attenuatton of vertical acceleration with depth, data 
normalized against acceloration at 5-foot depth, Nevada Test Site. 
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Figure 4.0 Attenuatlon of downward velocity with depth, 
Nevada Teal Site. 


ny 


sre, Sale 
ere 


ia 


a 


lags tie wave front by an increasing amount as the depth incrcases. This lag ts the sum of the 
overpressure rise thme plus a time lig due to spreading of the acccleration pulse with depth, 
For an Ideal overpressure wave form (no precursor), the latter time difference would corres- 
pond to the rise time of the vertic.:’ velocity pulse. The rise tlme (corresponding to an ideal 
overpressure wave form) bas been deduced from the Stanford Research Institute Priscilla ac- 
eeleration and velocity data (corrected for transmission time) and is shown in Figure 4.10. 
The single satisfactory obcervation on Tumbler ts also shown and compares surprisingly well 
with the Priscilla data. 

Still another method has bcen employed to show how velocity-time traces change with depth. 
The overpressure-time records and the velocity-time curves were normalized to their peak 
values; that is, cach paramcter magnitude was divided by the peak valuc and the ratio plotted 
versus time. Figures 4.11 through 4.15 present the results of these normalizations. Also in- 
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Figure 4.10 Increase in vertical velocity rise time with depth, Frenchman Flat. 


cluded in the figures are the normalized earth stress records, which will be discussed in Sec- 
tion 4.3.4. The records are aligned so that the peaks coincide and velocity curves are inverted; 
in this way easy visual comparisons can be made of rise times and decay times. 
Bue to the many overpressure gage cable breaks, Figure 4.11 indicates only limited conclu- 
sions. It appears that at Stations 1 and 2 (550-psi and 365-psi peak overpressure) the particle 
velocities at 5-foot depths attain a peak value in approximately the same time as the surface 
overpressure; however, at 10-foot depth the velocity rise time is about 5 msec longer. At 
Station 3 (Figure 4.12) the increased rise time for velocity is evident at both 5- and 10-foot 
depths and the decrease following the peak is close to that displayed by the surface overpressure. 
Similar behavior is observed at Station 4 (Figures 4.12 and 4.13) down to depths of 10 feet. But 
below 10 feet the velocity rise time increases significantly. Also the record appears to decay 
more slowly than does the surface overpressure. At Station 6 (Figures 4.14 and 4.15), where 
the precursor is well-developed, the rise to peak velocity fallows the surface overpressure 
down to 30-foot depth, but deviates marked!y at 50 feet. Also, at this station, one observes 
the first definite tendency for the velocity to decay more rapidly than the surface overpressure- = 
time measurement. At 2 depth of 50 feet (Figure 4.15}, the velocity rise time is about 20 msec 
longer than for the overpressure, which represents the only significant departure from follow- 
ing the input at Station 6; here again, however, the velocity decay Is more severe than the over- ° 
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Figure 4.11 Overpressure, particle 
velocity, earth stress, normalized 
to maximum value, Stations 1 through 
3, Shot Priscilla. 
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Figure 4.12 Overpressure, particle 
velocity, earth stress, normalized 
to maximum value, Stations 3 and 
4, Shot Priscilla. 
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Figure 4.13 Overpressure, particle velocity, earth stress, 
normalized to maximum value, Station 4, Shot Priscilla. 
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pressure decay. Station 7 (Figure 4.15) shows similar behavior to Station 6, there is only a 
small increase In velocity rise time at 10-foot depth. 


4.3.3 Displacement (Attenuation with Depth). The integrations of the Stanford Research 
Institute acceleration data are at present the only displacement data available for Priscilla. . 
These data are shown in Figure 4.16, normatized against the peak displacement at 5-foot depth. 
Additional data will be available for correlation when Sandla Corporation acceleration and long- 
ypan displacement gage data have been reduced (Reference 15). 

The attenuation of displacement for all ground ranges except the 100-psi station approximates 
the attenuation of velocity, Figure 4.9. At the 100-psi station the displacements below 5-foot 
depth are larger than would be expected from observations at the other stations. From the 
character of the acceleration records, ft is suspected that a refracted ground wave was super- 
imposed on the direct motion induced by the main air-blast wave. Due to the direct motion 
induced by the precursor, it is difficult to clarify this point since outrunning cannot be definitely 
established. Since outrunning occurs at depth prior to outrunning at the surface, ground wave 
interference ts rot unlikely and would explain the greater displacements, cf., Figure 4.33, 
Tumbler 1. The Project 1.5 displacement data expected on Priscilla may clarify this point. 

It should be noted that the displacement-time plots are probably not reliable beyond the time 
of peak value. The truth of this statement was demonstrated when, for the Project 3.5 data 
analysis (Reference 16) it was necessary to attempt correlation hetween the integrated dispiace- 
ments and scratch gages attached to a buried structure. Reading the acceleration-time records 
out to much later times and using the scratch gage data as a guide to proper corrections to be 
applied, it was found that the vertical displacements behaved differently from those shown in 
Figures 3.14 through 3.17. The primary difference was that, after reaching about the same peak 
value, the Project 3.5 integrations decreased steadily to almost zero displacement beyond 1.0 
second. The procedure used for Project 3.5 illustrates the distinct advantage of having an inde- 
pendent measurement with which to guide the corrections. 


4.3.4 Stress (Attenuation with Depth). The maximum vertical earth stress is plotted versus 
gage depth in Figure 4.17. For each gage station the curve has been joined to the maximum 
surface overpressure-datum point (plotted on the figure at 1-foot depth). Maximum stress, as 
determined by Carlson-type gages on Project 1.7, is also Plotted for the 1,050-foot ground range. 
There is some doubt about the validity of comparing Project 1.4 stress measurements with those 
from Project 1.7%. The latter gages were immersed in a large back-filled hole and the 20-foot 
deep gage was positioned on the bottom of the excavation. 

The attenuation of maximum vertical stress between the surface and 5 feet is not severe and 
is greatest at the first anr last station instrumented; the attenuation at these shallow depths is 
probably strongly a function of the duration of the overpressure input pulse. 

The marked increase in stress at 50 feet for both the 750- and 1,050-foot stations is unex- 
pected; however, there is little reason to expect these larger values of stress to be caused by 
over-registration of the gage. On the contrary there are many more situations which could 
cause under-registration of the Carlson stress gage than could cause serious over-registration. 

The pattern of attenuation of vertical stress with depth between 0 and 30 feet is in general 
agreement with the attenuation equation proposed on the basis of Upshot-Knothole data (Refer- 
ence 1), which predicts a halving of stress for each 10-foot increase in depth. Point by point, 
the data of Figure 4.17 depart severcly from this attenuation equation. However, the fit of data 
is only slightly poorer than the original data on which the equation is based. The 50-foot 
Stresses, however, tend to deny the validity of this attenuation equation. 

Referring to Project 1.4 and Project 1.7 results plotted on Figure 4.17, the large discrepancy 
between the maximum Stress at 5 and 10 feet, measured at the two adjacent locations, is readily 
apparent. Since no systematic behavior is observed at the two locations, it is difficult to reason 
that the differences are due only to the different plantin:: procedures used by the two projects 
(also see wave forms of Figure 3.23). Since these deviations represent wide variations in stress 
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Figure 4.14 Overpressure, particle velocity, carth stress, 
nornmulized to maximum value, Stations 5 and G, Shot Priscilla. 
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Figure 4.15 Overpressure, particle velocity, eurth stress, 
normalized to maximum value, Stations 6 and 7, Shot Priscilla. 
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Figure 4,16 Attenuation of vertical displacement, Frenchman Flat. 
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Figure 4,17 Maximum earth stress versus gage depth, Shot Priscilla. 


response, one ts not very encouraged about the uscfulocss of a limited number of Carison gage 
results in predicting ground motion and/or structure response. 

One cause of decay of peak stress with depth Ls the finite duration of the blast wave on Shot 
Priscilia. An inspection of the air-blast wave forms shows that the duration of the main peak 
of the blast wave was relatively short. Measuring between points of arrival and one helf the 
peak pressure, this was about 28 msec at 750 fect and 55 msec at 1,050 f[cet. By rough calicu- 
lations of trace velocity, these correspond to a wave length of half-peak pressure of only about 
100 to 125 feet. This ts certainly not long cnough to be considered infinite compared with gage 
depths of 50 or even of 20 fect. There is reason to believe that, had the duration been much 
longer, the stresses measured at 20 and 30 [eet would have more nearly approached the pcak 
air-blcst pressure. 

There are several ways of looking at the effect of changing soil characteristics upon the at- 
tenuation with depth of maximum vertical earth stress. 

One is to assume initially that the gages record true {ree-field phenomena and then look at 
the effect of the hard and soft layers upon wave reflection patterns. By so doing, one can give 
a possible explanation for an increase (i.e., negative attenuation) in response with depth. How- 
ever, this argument applies only if the duration of the peak overpressure is short compared 
with the time required for the wave to be reflected once again at the earth’s surface. Moreover. 
it is Lmprohahle that the high measured stresses recorded at 50 feet could be explained on this 
basis; the stress-time curves are relatively flat in the region of maximum stress and give no 
indicatior of reflections. Furthermore, the Project 1.7 20-foot gage was at the bottom of the 
backfill and did not exnibit an increase in stress. 

Another way to approach the problem is to consider it statically and to surmise how the load 
at the surface was carried down through the soll. During the backfilling of the deep holes at 
Stations 4 and 6, a geophone was planted in each hole near the wall at a depth of 50 feet, and 
measurements were made of the travel time from the surface to each geophone. These travel 
times correspond to average velocities from 0 to 50 fect of 1,390 + 40 ft/sec and 1,450 + 60 
ft/sec, compared with 2,060 ft/sec from Shot Point» 1, 2, and 3, (Figure 2.6). These data 
tmply that the backfill failed to reach the original modulus and that the effect of this condition 
reached somewhat beyond the edge of the hale. 

Qn the assumption that the stress gages were placed in a medium of uniformly low elastic 
modulus, one would expect the measured stresses to be lower than those of the free field, with 
the stress differences largest where differences in elastic modulus (or seismic velocity) are 
largest. This conclusion is based on the premise thar the boundary between lower and higher 
elastic modulus media can carry vertical shear. This assumption is good for the undisturbed 
medium and probably satisfactory even at the backfill-hole wall boundary. Under loading, the 
two media will undergo much the same strain; hence, the stress will be lower in the media of 
lesser modulus. 

Note, however, that the local seismic velocity at 40- to 50-foot depths in the free field (Fig- 
ures 3.25 and 3.27) is only slightly higher than that which appears characteristic of the average 
at the 750- and 1,050-foot stations. This fact may offer a possible explanation, in the light of 
the above discussion, of the higher measured stresses at both 50-foot depths. RB can be con- 
eluded that the 50-foot stresses are more indicative of th: free-field stresses and that those 
measured at intermediate depths were depressed due to the lowered modulus. 

Still a third way to view the problem of explaining the stress results is to take a somewhat 
eyaical approach. One can say that the backfill in the niche around the stress gages at 20- and 
30-foor depths was not well compressed, anc hence these gages saw only a minute stress. 

To summarize, it is probabic that all these effects are present in varying degrees of severity. 
Sullice it to say chat it is difficult to dctermine the validity of the {ree-field carth stress results 
obtained by Project 1.4 0n Priscilla. The control of backfill operations during Priscilla was 
considered optimum and extraordinary care was taken in placement of Carison gages; still the 
Stress measurements appear inconsistent internally and hence leave much to be desircd. 2 is 
believed af Ubis time that the present concept and method of earth stress measurement may be 
inadequate and misicading. Alternative methods cannot as yct be offered, except that probably 
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stress must be derived from other primary measurements such as strain and acceleration, 
which are more Independent of hole size and the character of the backfIlL 

Reference to the normalized records of Figures 4.11 through 4.15, reveals some intcresting 
facts about the stress-time wave form changes with depth. At Station 3 (Figure 4.12), the rise 
time to peak stress increascs noticcably from §- to 10-foot depths; hnwever, the decay of stress 
at both depths, ts comparable with the decay of the overpressure input. Figures 4.12 and 4.13 
include the results at Station 4; here no significant increase in streys rise timc is noted at 1-, 
$-, and 10-foot depths. The stress records at 20- and 50-foot depths deviate markedly from 
the input overpressurc-time wave form; also, they do not resembic cach other. 

At Station 6 (Figures 4.14 and 4.15), the stress rive time and decay follows the overpressurc 
input quite closely to depths of 10 feet. At 20 fcet, the rise time is noticcably longer for the 
stress, whereas at 30 fect the stress record Is strange, and at 50 fect the stress response to 
the precursor portion of the Input wave is depressed relative to the peak. The one comparison 
possible at Station 7 reveals that at 5-foot drenth the stress record followed the overpressure 
inpat quite faithfully. 


4.3.5 Strain (Attenuation with Depth). Due to the full or partia! failure of eight of the ten 
short-span strain gages, the strain-time data of Project 1.4 were rather mozger. On some of 
the gartial records it was possible to Identify the maximum vertical strain recorded at the gage 
station; thus, the measurements ytelded a total of five peak strain values, three at 750-foot 
ground range and two at 1,050 feet. 

ln addition to the strain gage measurements, it was possible to compute the average vertical 
strain from the integrated accelerations (L.c., displacemert-time plots). Using the strain-time 
results as 2 guide to give the time of maximum strain; these computations were carried out and 
are listed in Table 4.2. The results are plotted in Figure 4.14 and show that the calculated max- 
imum average strains are not inconsistent with the few strain-time measurements obtained. 
This takes the form of an independent check and leads to more confidence in the manner in which 
the acceleration-time records were integrated. 

The data of Figure 4.16 indicate that at the 750-foot station (275 psi maximum overpressure) 
the peak earth strain bad 2 decided atternzation with depth, with a pronounced decrease between 
1 and 30 feet below the ground surface; at larger depths (to 50 feet), the maximum vertical strain 
appeared .» level off at about 5 parts per thousand. On the other hand, the data obtained at 1,050 
feet (100-psi overpressure level) showed practically no change in maximum vertical strain with 
increasing depth; In fact, there was even a slight indication of larger (negative attenuation) strain 
near the 30-foot depth. The reason for this difference in behavior was probably traceable to the 
character of the input overpressure-time 2t the two stations; this is discussed in more detail in 
Section 4.4.4. - 

Tt is noted that the Station 4 strain measurement at 5-font depth (4SV5A) had a much smaller 
peak magnitude than indicated by the computations using integrated accelerations. Also, the 5- 
foot measurement was even smaller at maxima than was the 10-foot strain at the same station 
(4SV10). Another pertinent point can be made upon inspection of the 4SV5 gage record (Figure 
3.24). The record indicates that although the overpresscre input increased from a precursor 
pressure of about 25 psi to a peak pressure exceeding 225 psi, a nine-fold increase, the strain 
merely doubled. The 4SV50A and 6SVS0 strain records are examples of characteristic behavior 
in this regard. 

tt must be concluded, therefore, that the 4SV5 strain measurement was, for some reason, 
erroneous; the true peak strain at this position was significantly larger than was recorded. 


4.4 GROUND SHOCK AND OVERPRESSURE 
4.4.1 Acceleration (and Overpressure). To compare accelcrations on the basis of air-blast 
input, the ratlo of pe2k accclcration to peak overpressure Versus overpressure was plotted. 
Figure 4.1$ summarizes all 5-fvot depth vertical acccleration observations in Frenchman 
&2 
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Flat. Also plotted are the accclcration-overpressure ratios for the Priscilla and Upshot- 
Knothole Shot 10 precursor. The measured overpressures at Priscilla Stations 2 and 4 2zgeared 
low when compared with pressure dat2 of previous shots. At Statlon 2, a cable break occurred 
prior to recording of the peak overpressure. The source of the apparently low pressure at Sta- 
tion 4 has not been found. Both the as-read and smoothed data points are indicated. The accel- 
erations induced by the main alr-blast wave decayed rapidly on Priscilla due to the formation of 
the precursor and also due to the distance effect (Tumbler data discussed below). Taking this 
into account, the Upshot-Knothole Shot 10 Type 1 overpressure wave form accclerations were 

in agreement with those of Priscilla. The single Upshot-Knothole Shot 9 observation appeared 


TABLE 14.2 MAXIMUM STRAIN, SHOT PRISCILLA 


Computations Measurements 

Gage . Disp] Aver Max Gage Peak Time 
Code Perth Times ts ort. Strain Code Strain of Peak 

It sec ft ppk* pek® sec 
4V1 1 0.193 0.390 23.8 — — _— 
4vs 5 _ 0.295 _ — a oe 
4vi 1 0.195 0.411 26.6 _— — — 
4vVilo 10 _ 0.171 — 4SVS 8.8 0.195 
4V5 5° 0.205 0.409 30.0 _ — — 
4Vlo 10 —_ 0.259 _ 4SV10 13.7 0.215 
4Vi0 10 0.225 0.365 1L2 _ _ _ 
4V20 20 -_ 0.253 —_ _— — — 
4V20 20 0235 0304 7.0 ~— —_— — 
4V30 30 _ 0.234 _ —_ — — 
4V30 30 0.220 0.254 63 _ — —_ 
4vso 50 _ 0.127 _ 4SVSa 7B 0.249 
evi I 0.284 0.147 5.8 _ — _ 
6VS 5 — 0.124 ~ _ —_ _ 
6VvS 5 0.250 0.156 41 _ _ _— 
6V20 20 _— 0.080 _ _ — —— 
6V20 20 0.300 0.122 6.5 _ — — 
6V30 30 _ 0.0352 = 6SV20 59 0.304 
6V30 30 0.317 0.110 25 _— — — 
6V50 50 _ 0.061 —_ @svso 20 O7390 


© Perts per thousand. 


lower than would be expected from Tumbler Shot 1, and the Priscilia precersor accelerations 
were larger by about the same amount. In general, the accelerations appeared to be larger in 
the Frenchman Flat area than in the Yucca area, not however by an amount which was significant 
compared with the scatter of data. 

Cf some interest were the peak acceleration-overpressure ratios obtained [rom the horizontal 
(positive outward) acceleration observations at 10- and 50-foot depths. At Station 4, the ratios 
were about one third those corresponding to the vertical accelerations and at Station 6, although 
the 10-foot depth vertical gage was lost, the 50-foot data indicated that the horizontal and vertical 
components were approxzimatcly cqual. 2 

The ratio of peak negative acceleration to peak overpressure Versus overpressure is plotted 
in Figure 4.20 for Tumbler Shots 1 and 2 (References 2, 17 and 18). The ratio appears to be 
fairly independent of overpressure. The heights of burst of both Shot 1 and Shot 2 were sufli- 
ciently high that no precursor formed and all overpressure wave forms were nearly iccal. 
Since the yields were nearly identical, it must be concluded that the two to one difference in 
acceleraiinn-pressure ratio must result from differences in soil properties between the French- 
man Flat and Yuce2 test areas. 

Figure 4.21 compares the three Tumbler shots detonated over the same portion of the Yucca 
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test areca. One first notes the marked reduction in the Shot 4 acceleration-pressgure ratio at 10 
psi. The smuiler than expected accelerations were the result of the nonideal precursor wave 
forms. The Type 4 overpressure wave form has the largest rise time and hexce induces tke 
least ground acceleration. As soon as the wave forms return to near ideal (Type 7}, the accel- 
eration assumes its normal value. Taking the effect of overpressure wave form into considera- 
tion, the effect of device yield if any, dues not appear to be discernible within the scatter of data, 
In Figures 4.20 and 4.21 there is an obvious tendency for the acceleration-pressure ratio to 
increase at both ends of the pressure scale. When the data of Figure 4.21 are plotted against 
ground range as In Figure 4.22, there appears to be a fair correlation with distance. This 
characteristic may be due to changes in soll properties or due to shock wave parameters not 
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Figure 4.18 Earth strain versus depth, Shot Priscilla. 


considered, e.g., shock velocity and/or wave front curvature. The increase in acceleration 
ratio at close-in stations may poss-bly be connected with the latter since the same tendency is 
observed on Tumbler Shot 1 (Figure 4.20) and on the two near ideal wave forms on Priscilla 
(Figure 4.19). This point cannot be clarified at the prescat time due to insufficient theorctical 
investigation and without more detailed study of soil propertics in the Yucca arca. 

The remaining sct of observations at 5-foot depth are the vertical acceleration data of Jangle 
HE-4 (Referenec 4), 2,560 pounds of TNT detonated at the ground surface. The acccleration- 
overpressurc ratios are plotted in Figure 4.23 and are consistent with previous observationc. 

Jangle (Ref-:rence 19) vertical acccleration data are shown in Figure 4.24. The scatter of 
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Figure 4.19 Vertical acceleration, Frenchman Fiat, $-foot depth 
(numbers indicate overpressure wave form). 
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Figure 4.20 Vertical acceleration, Tumbler Shots 1 and 2, 5-foot depth 
(arrow points to first ground wave outrunning wave form). 
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Figure 4.21 Vertical acceleration, Tumbler Shots 2, 3, and 4, 5-foot depth 
(numbers adjacent to data indicate overpressure wave form). 
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Figure 4.22 Vertical acceleration at 5-foot depth versus ground range, Tumbler 
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these data ls much greater than on any of the previous experiments discussed, which makes 
analysis of these data rather tenuous. The attenuation indicated by the Jangle S 10-foot data 
appears to be greater than that observed on Tumbler or Priscilla, as was discussed carller. 

The vertical acccleration data of Project Mole (Reference 3) aboveground detomtlons are 
shown in Figures 4.25 and 4.26. Nevada sand and gravel mix corresponds to the soll tn the 
Jangle S, HE-4 test area. The Utah dry clay site was located at Dugway Proving Ground, Utah. 

In Figure 4.25, although the scatter of data is Large, {note that acceleration-pressure ratios 
are plotted on logarithmic scales), no systematic variation with height of burst is to be found. 
Secondly, the range over which the acceleration-overpressure ratio extends is an order of mag- 
nitude large, i.e., the ratio is more dependent on overpressure (or ground range) than that of 
the large HE (Jangle HE-4) or nuclear (Jangle S) detonations in the same area. This trend to- 
ward higher ratios at the extreme ground ranges is observed however in Tumbler Shots 2, 3, 
and 4 and may be connected with outrunning of the ground wave. 

tn Figure 4.26 the dependence on overpressure is not as pronounced, and the accelerations 
appear to be a factor of two larger than those of Figure 4.25 for pressures greater than 30 psi. 


4.4.2 Velocity (and Overpressure). Subsequent to outrunning of the ground wave, the velocity- 
jump overpressure ratio increases with decreasing overpressure (cr increasing ground range). 
This observation is in conformity with elastic theory (private communication with J. K. Wright, 
AWRE). Prior to outrunning, the ratio should be independent of pressure. All 5-foot depth 
peak velocity data, where the air blast was superseismic, are shown in Figure 4.27; 80 percent 
of the data fall within a velocity-pressure ratio range of 0.04 to 0.06. There appears to be no 
systematic variation with either yield, overpressure level or wave form, or test area. The last 
conclusion is at variance with the behavior of the acceleration data. 

The data shown in Figure 4.27 may be compared with the theoretical result, 


Uy 1 _ 32.2 x 144 = 0.032 ft/sec 
P "pc, © 90x1,600 ~ " Psi 


Where the soil density is taken at 90 pcf, 2nd from Figure 3.30, Cz, at 5-foot depth is 1,600 
ft/sec. The agreement between theory and experiment appears extremely good, especially if 
one recognizes the approximate nature of the theoretical result. Secondly, the true seismic 
velocity in the first 5 feet is difficult to measure accurately: the estimated error is of the order 
of 30 percent. The Large seismic velocity gradient also presents both a theoretical and practi- 
eal problem. 

The velocity of propegation of ground motion may be determined directly during an experi- 
ment from the arrival time of ground motion and the air blast arrival over the gage location. 
This calculation was made for Priscilla using arrival times at 10-foot depth. These calculations, 
Table 4.3, lead to the conclusion that the average propagation velocity of ground motion between 
the surface and 10-{foot depth is less than 1,000 ft/sec. HM the average velocity of Table 4.3 is 
used in the previous calculation one arrives at the result, 


Uy 
= 0.060 

The fact that the high amplitude stress wave propagates slower than the seismic wave is not 
unexpected in view of the theoretical treatment of von Karman and Duwez (Reference 20) of wave 
propagation in plastic solids. However, additional treatment of this observation is required be- 
fore definite conclusions can be derived. 

Velocity-jump-overpressure ratios for Tumblcr Shots 1 through 4 are shown in Figures 4.28 
through 4.31. All data displayed tn Figure 4.28 have been reintezrated as outlined in Appendix 
C. Data in Figures 4.25 through 4.31 have been reintegrated only for those stations where the 
blast wave is superscismic (the first station at which the ground wave outruns is also included). 
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Figure 4.23 Vertical acceleration, Jangle HE-4, 5-foot depth 
(f indicates interpolated overpressure data). 
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Figure 4.24 Vertical acceleration, Jangle S. 
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Figure 4.25 Vertical acceleration, Project Mule, 
Nevada sand and gravel mix, 5-foot depth. 
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Figure 4.26 Vertical acceleration, Project Mole, Utah dry clay, 5-foot depth. 
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Figure 4.28 Vertical velocity, Tumbler Shot 1, 5-foot depth. 
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The remaining data arc the original reported results. In Flgure 4.32 are shown all 5-foot ver- 
tleal velocity observations tn Frenchman Flat. 

The horizontal outward velocity at Station 4 on Skot Priscilla gave a peak value at 10- and 
50-foot depths which !s about onc sixth the maximum downward velocity at the same locations. 
Although the 10-foot vertical accelerometer record at Station 6 was lost, the comparison at 50 
fect reveals that the maximum outward velocity was less than one half the peak downward vc- 
locity. The apparent relative increase in the outward component a3 one goes to increasing 


depths and ground ranges is due to the greater Influence of signals from remote sources upon 
the ground motion, 


4.4.3 Displacement (and Overpressure). The peak vertical displacement-positive overpres- 
sure impulse ratios at 5-foot depth for all records analyzed to date are shown in Figure 4.33. 
Some remarks on the Priscilla results are periinent. First, the 5-foot depth displacement at 


TABLE 4.3 AVERAGE PROPAGATION VE- 
LOCITY OF DIRECT WAVE, 0 
TO 10 FEET, SNOT PRISCILLA 


Overpressure cL 

Fi ft/sec 
550 685 
+40 660 
340 3980 
278. 1,080 
210 ISS 
100 880 

60 39380 


Average value 960 


Station 2 has been estimated from the displacement at 10-foot depth, Point e of Figure 4.33, 
and the attenuation of displacement with depth (Figure 4.11). Second, the overpressure and 
impulse at Station 2 have been estimated by a reasonable interpolation between Stations 1 and 3. 

The smaller-than-expected displacement at Station 1 on Priscilla has been checked and is 
further verified by the 10-foot depth displacement. 

Considerable variation exists between the Tumbler displacement impulse ratios. No further 
correlation appears evident through consideration of device yield since the variation of displace- 
ment impulse ratio is not systematic. The agreement between Tumbler 4 and Priscilla is en- 
couraging but possibly fortuitous. 

Figure 4.34 displays reed-gage displacement data for a number of Plumbbob detonations in 
the Yucca test area (Reference 21). Positive impulse values have been estimated from meas- 
ured overpressures and device yield using the correlation presented in Appendix D. The ver- 
tical lines adjacent to each datum point indicate the passible error involved in this estimate. 
Corresponding displacements have been computed from Priscilla accelerograms (Section 4.5). 
The same degree of variation is evident between these data as was found in Figure 4.33. 

The reasons for the wide scatter of displacement data are not known. The only conclusions 
which appear valid at the present time are that approximately the same degree of variation ex- 
ists in Frenchman Flat as between the Frenchman Flat and Yucca test areas, and that displace- 
ments twice those observed on Priscitla are possible for similar geological conditions. 

The variation of displacement-impulse ratio with peak overpressure is not surprising phe- 
nomenologically; however, fram an empirical prediction standpoint it would be desirable to 
climinate such a dependency. Several inquirics have been made along these lines, and it ap- 
pears possible that displacemert may be proportional to that fraction of overpressure impulse 
above a certain pressure or to that fraction of overpressure impulse delivered prior to a cer- 
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Figure 4.29 Vertical velocity, Tumbler Shot 2, 5~foot depth 
(f indicates as-reported data). 
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Figure 4.30 Vertical velocity, Tumbler Shot 3, 5-foot depth 
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Figure 4.32 Vertical velocity, Frenchman Flat, 5-foot depth 
(numbers adjacent to data indicate overpressure wave form). 
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Figure 4.33 Summary of vertical displacement, Nevada Test Site 
(e indicates estimated value). 
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tain fixed time duration. Additional work in this arca is requircd before any conclusions can 
be presented, 

In any analysis of vertical displacement, eximination of the validity of the double integration 
of acceleration is worthwhile clfort. As previously mentioned, most first integrations of accel- 
eration lead to nonzero velocities at the end of the record. Provided the velocity baseline shift 
is linear and onc ts able to make a rellable estimate of the time of zero velocity, the second 
integration affords little difficulty. Prior to Priscilla, however, data did not exist which would 
serve as an independent check on the results obtained. The Priscilla long-span strain gage data 
now provide this opportunity. In Table 4.4 a comparison is made between the Sandia long-span 
strain gage results and the Stanford Research Institute double integrations. Only near-surface 
data are compared since strain gage data at depth have not yet been reduced. The agreement 
between these two Independent measurements is surprisingly good, except at 1,350-foot range. 
One can speculate that the Sandia anchor at 200 feet at this station would be most Lafluenced by 
remote source signals, thereby reducing the apparent near-surface peak displacement. How- 
ever, this type of reasoning is dangerous without more complete data. It should be noted, how- 
ever, that the capability of the integration method to lead to valid permanent displacements is 
severely limited. 


4.4.4 Strain (and Overpressure). Reference to Figure 4.18 indicates quite a differcnt peak 
strain attenuation with depth at the two stations instrumented (750- and 1,050-foot ground range). 
The reason for this difference may be traced to the different character of the inputs at the two 


TABLE 4.4 COMPARISON OF MAXIMUM 
TRANSIENT-VERTICAL DIS- 
PLACEMENT, SHOT PRISCILLA 


Maximum Displacement 
Ground Range = Spi, 5 ft* Sandia, 3 ft? 


t fc ft 
650 0.84 0.82 
850 >0.405 9.52 

1,050 0.27 0.23 
1,350 0.24 Q.12 


* Double integration of acceleration. 
+ Long-span displacement gage. 
} Cable break, displacement estimated. 


Stations. At Station 4 (750-foot ground range), the precursor pressure of about 25 osi precedes 
a rise (in about 5 msec) to 228-psi peak overpressure. By contrast, the input at Station 6 indi- 
cates a precursor pressure of 20 psi and rises to only 104-psi peak in about 15 msec. 

Probably the rise time associated with the input overpressure assumes the greatest impor- 
tance when considering earth strain, particularly near the ground surface. Measurements of 
Seismic velocity indicate that the weathered layer at Frenchman Flat extended to depths of 5S or 
more feet. It is feasible that an overpressure-time input such as 4B (Figure 3.1) would pro- 
duce large strains in such a layer, whereas the 6B input would not. The few strain measure- 
ments obtained by Project 1.4 severely limit firm conclusions concerning the phenomenology. 


4.5° RESPONSE SPECTRUM OF GROUND MOTION 


4.5.1 Theory. The response spectrum is defined as the maximum response of a linear, 
single degree of freedom, Spring mass system, relative to the motion of ground. Sucha sys- 
tem is shown in Figure 4.35. 

The equation of motion, in terms of the relative displacement, x, {ts (dot superscript denotes 
differentiation with respect to time, } 


Ke Onwk + wrx ~alr) 
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Where: n Is the ratio of damping to critical damping 


@ = undamped natural frequency, w = k/m 
m = mass of system 

k = spring rate of system 

a = ground acceleration 

t = time 


Development of this equation may be found in any standard text on dynamic motion, or in Ref- 
erences 22 and 23. 


Let y = x and integrate once, resulting in 
¥ + 2nwy + wy =-u(t) 
Where u ts ground particle velocity. Let z = y and integrate once again, 
Z+mwz—-atz =-s(t) 


Where s is ground displacement. 


These equations demonstrate that the response may be calculated from knowledge of any one 
of the ground motion parameters, acceleration, velocity, or displacement. 
Three response spectra may be defined: 


Ixl , the displacement (response) spectrum 
max 
witlmax, the velocity (response) spectrum 
wixlmax, the acceleration (response) spectrum 


The velocity spectrum is proportional to the maximum strain energy per unit mass for an 
instantaneous applied force and no damping, 


2 klix! 1 2 
& . —=Mmaz .«£—- 
=—=-3 3 wx | ax 


For a smali amount of damping, n~ 0, the acceleration spectrum equals the maximum 
absolute acceleration of the mass, q. 
One has 


x = q-sorx = q-a 
From the equation of motion(n = 0) 
X+2 =-w"x. 
Hence 


lalmax = @Ixhnax- 


A few characteristics of response spectrum may now be lUustrated. Take the second form 
of the equation of motion. The ground particle velocity may be characterized by 2 peak value, 
Um, and a duration T, as In Figure 4.36, i.e., u/Ug, = (7) andt = T>. 


Let ¥ = (y/u,,)T* and 2 = wT, then Y + NY + N?Y =-&7). 
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The velocity spectrum, &lY¥l.4,, 1s a single valued function of 2 only. 
But -* 


: 1 _ ixlmax 
aly. = wT oe er ae ® 


Hence (wlxl..45)/tm is also a single vaiued function of wT as in Figure 4.372. 

In Figure 4.37b the velocity spectrum of Figure 4.36 is displayed for large T, large yield 
devices and small t, small yield devices. The peak value of wlxl,..,, does not change pro- 
vided the peak particle velocity does not vary, i.e., constant overpressure. However, the 
curve shifts toward the lower {rcquencies as the yield increases. These statements are true 
only if $(7) is invariant. 

lh the same manner {t may be shown that 


Ixl max 


an = f; (wT) 


wixlnax = y, (wr) 


2m 


Where d,, and am are the peak ground displacement and acceleration, respectively. 

As the frequency w (or more rigorously wT) approaches zcro, the response |xl,.,, approaches 
the peak ground displacement. At large values of w (wT) the response approaches the peak accel- 
eration divided by u*. Hence 2s w— > 0 


1 
Ixlmax~dm = Ym T — dr 
0 


pri 
and 
wlx!max. C, eT 
Um 
where 
1 
Cc, = ar 
6 Um 
Asw—>o 
a/(u/um) 
atx! ~a, = Sm jun! 
max ~ “m "pT Vas | max 
and 
wlxl max c, 
Um wT 
where 
CG = d/(u/um) 
° dz max 


The conditions undcr which the velocity spectrum may be scaled as in Figure 4.37b, require 
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Figure 4.35 A linear, single-degree-of-freedom system. 
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Figure 4.36 Schematic diagram of ground particle velocity. 
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Figure 4.37 Schematic diagram of velocity spectrum. 
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that the maximum displacement at ccnstaat overpressure (u,, |S proportional to overpressure) 
vary as T and that the peak acccleration at constant overpressure vary as T™'. That is, all 
details of the velocity time history must scale tn 2 consistent manner. 


4.5.2 Measurements. The displacement spectrum may be determined cither by direct meas- 
urement with reed gages, by numerical calculation of the response of a lincar system to the 
ground motion, or through use of the electrical-mechanical analogy by converting the measurcd 
ground motion into an electrical signal (Reference 24). All methods will produce identical re- 
sults provided certain precautions are taken. An excellent example of equivalent results ob- 
tained by direct measurement and the analog calculation is presented in Reference 25. 

Reed gages should be constructed so that their springs remain linear up to the maximum 
response measured. The damping should be linear, ie., viscous, up to the maximum response. 
The percent of critical damping should be known and preferably should be small and approximately 
the same for each spring of the set. 

The accuracy of numerical or analog calculations of the response depends first on the accu- 
racy of measurement of the ground motion, i.e., the accelerogram. As discussed in Appendix 
C an acceierometer-galvanometer system will yield faithful results if its natural frequency is 
more than twice the fundamental frequency observed in the ground motion, and hence it is un- 
wise to calculate response spectrum for frequencies greater than half the accelerometer or 
galvanometer frequency, whichever is smaller. — 

Secondly, the accuracy of numericat calculation depends on the fineness of the calculation 
network, Le., the time interval between successive readings of the record, in relation to the 
sophistication of the integration program. The time interval is usually limited by the precision 
of the time base on the original record (0.5 msec on Priscilla, 1 msec on previous operations). 
The integration program selected is a compromise between the cost of running a complex com- 
putation and the maximum [frequency to which the computations are carried. For example, the 
Simplest program, using the particle velocily as Une input, appears to vive reliable results for 
the majority of the Priscilla records to a frequency of 250 radians/sec (40 cps). To compute 
the response of a 100-cps system, a more complex (and costly) program must be written. 

This second consideration is not usually a limiting factor in analog calculations. However 
a necessary condition for both calculations is that the velocity at the end of the pulse be zero. 
This is of major importance at low frequencies and becomes less important 2s the frequency 
increases (Reference 26). The importance of correcting for the baseline shift discussed in 
Appendix C is now further emphasized. 

Response spectra for all the Priscilla ground motion data, excepting Stations 2V5, 5V5, and 
SV10 have been calculated to a frequency of 150 radians/see (24 cps) and about half these have 
been extended to 250 radians/sec (40 cps). Station 2V5 was eliminated due to a spurious elec- 
trical signal, and 5V5 and 5V10 were eliminated hecause of cable breaks. Displacement response 
specbra are shown in Figures 4.38 and 4.39 and the corresponding velocity response spectra in 
Figures 4.40 and 4.41. In Figures 4.40 and 4.41 are shown the asymptotic values as w--->0 
(Curve A) and as w---><: (Curve B). In Figure 4.42 are shown a selected number of accelera- 
tion spectra. The calculations have obviously not been carried to frequencies high enough to 
determine the maximum of the acceleration Spectra. 

The results of Figure 4.38 indicate a flat response to 10 radians/sec (1.6 cps). A reed gage 
frequency of 2.5 eps would respond generally within 10 percent error in maximum transient 
vertical and horizontal displacements. The spectra are smooth out to and including Station 4 
(275 psi). The spectra at Station 6 (100 psi) begin to assume a jagged character which becomes 
more pronounced at Station 7 (60 psi) as the frequency increases, 

The overpressure wave farms from Statior 3 (340 psi) on were precursor Type 1, Appendix 
A. The only significant difference in the shape of the overpressure-time curve was increased 
precursor duration and rise time of the main overpressure pulse at the larger ground ranges. 
Herce, it does not appear that the change in character or the response can be associated with 
the local ground wave but rather that it results from interference of the refracted ground- 
transmilted waye. 
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Figure 4.38 Continued. 
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Figure 4.39 Shot Priscilla horizontal displacement spectra, 
0.5 percent critical damping. 
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This conclusion is supported by close eximinitlon of the acceleration records. At Station 4 
(275 psi), the ground wave has not yct outrun the main air-blast wave but arrives slightly behind 
the direct signal, The interference ts most pronounced at the 20- and 30-{vot depths. At Sta- 
tions 6 (100 psi) and 7 (60 psi), the ground wave appears to be outrunning the main blast wave 
but not the precursor. This is possible since the precursar propagation velocity at these ground 
ranges is roughly twice the min wave velocity. 

This example illustrates one of the inherent advantages of accelerograms over recd-gage 
data in that significant information may be lost when an attempt is made to reduce primary data 
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Figure 4.40a Vertical velocity spectra, 5~foot depth, Shot Priscilla. 


at the instrument itself. However, the recd gage, or preferably a displacement-time meter, 
has the advantage of being insensitive to the type of baseline correction which appears necessary 
to the integration of accelerograms. These comments point out that neither accelerometers nor 
reed gages (or other displacement meters) are completely self sufficient and that the most ad- 
vantageous instrumentation probably will consist of both types placed as a pair rather than as 
Separate entitles. 


4.5.3 Correlation and Scaling. The question of correlation and scaling of the response spectra 
is now considered. In Figures 4.43 and 4.44 are shown the vertical velocity spectra for 5- and 
10-foot depth normalized in the manner indicated in the section on theory. Data fcr the 100~- and 
60-psi stations arc not shown duc to the interference of the ground wave, and data for the 550-psi 
station are not shown due to the anomalous (low) values of peak displacement at this station. 

The velocity spectra were chosen over the displacement and acceleration spectra since the 
wlximax, i.c., the maximum energy per unit mass, Ls one of the principal factors in design of 
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Figure 4.40b Vertical velocity spectra, 16-foot depth, Shot Priscilla. 
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Figure 4.40c Vertical velocity spectra, 275 psi, Shot Priscilla. 
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Figure 4.40d Vertical velocity spectra, 100 psi, Shot Priscilla. 


shock tsolation equipment. Since the displacement and acceleration are both known functions of 
the velocity, they may be easily derived from Figures 4.43 and 4.44. For T, the normalizing 
factor for w (Section 4.5.1), the overpressure positive phase duration is used, again as a matter 
of convenience, rather Lhan the duration of the velocity pulse. Whether or not this serves to 
correlate the data consittutesy a test of the validity of the procedure. Furthermore, the over- 
pressure duration ls a known quantity independent of the ground motion, a factor necessary for 
sealing to larger yields. 

The smoothed curves of Figures 4.43 and 4.44 are identical, which is encouraging; however, 
the variation in T ts only 25 percent and hardly constitutes a valid test. 

The normalized velocity spectra for 50-foot depths are shown in Figure 4.45 and compared 


Ix]ag. — inches 


X\X 


Figure 4.41a Horizontal velocity spectra, 275 psi, Shot Priscilla. 


with the spectra for 5- and 10-foot depths. The normalized value does not change significantly 
although the frequency at which this maximum occurs decreases. Table 4.5 summarizes the 
maximum ratios for all gages. The average ratio for all vertical gages except Station 7 (60 psi) 
is approximately 1.2; at Station 7 the ratio Is greater than 2 due to the ground wave. Note that 
for all horizontal gages the maximum value of the normalized velocity is greater than two. 

The normalized horizontal spectra for the 10- and 50-foot depths are shown in Figures 4.46 
and 4.47. The 100- and 275-psi, 50-foot depth and the 100-psi, 10-foot depth spectra are sim- 
ilar while the 275-psi, 10-foot depth exhibits marked differences. 
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Figure 4.42 Shot Prisoiila vertical acceleration spectra, 
0.5 percent critical damping. 
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Figure 4.43 Normalized vertical velocity spectra, Priscilla 
5-foot depth, 0.5 percent critical damping. 
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Figure 4.44 Normalized vertical velocity spectra, Priscilla 
10-foot depth, 0.5 percent critical damping. 
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Figure 4.46 Normalized horizontal velocity spectra, Priscilla 


100 
10-foot depth, 0.5 percent critical damping. 
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4.6 SOIL STRESS-STRAIN CONSIDERATIONS 


4.6.1 Soil Survey Tests. As a part of Project 3.8 carried cut by Waterways Experiment 
Station (Reference 14), the compressibillty characteristics of the undisturbed foundation were 
determined using [our types of tests. A modulus of deformation was determined from consoll- 
dation test data using the slope of the stress-strain curve in the first cycle of loading at applicd 
Stresses of 50 psi and 100 psi. The modulus of compression was determined using the maximum 
Slope of the stress-strain curve [rom constant ratio of applicd stress triaxial tests. A compres- 
Sive modulus was determined from soniscope test data using conventional procedures. A modu- 


TABLE 4.5 MAXISIU31 OF NORMALIZED VELOCITY 
SPECTRUM, SHOT PRISCILLA 


wlx} > Umax 
Depth Ovecpressure a 


Vertical Horizontal 
ft psi 

I 275 1.32 _ 
200 1.23 —_ 
$s 550 1.22 _ 
1.34 —_ 
275 1.22 _ 
100 1.50 _ 
6 0.98 — 
10 350 L30 _ 
440 1.18 _ 
wo 13 _— 
235 1.08 2.0 
is] _ 2.0 
E ao 1.52 _ 
20 203 L3 _ 
100 1.31 —_ 
30 27s 1.46 _ 
100 1.18 _— 
50 275 1.17 1.9 
- 14g 1.28 23 


lus of elasticity was determined using data from field plate bearing tests. The data from the 
above tests which are pertinent to Project 1.4 are summarized in Table 4.6. 

Looking at Table 4.6, based upon the consolidation test, it may be noted that the modulus of 
deformation for the natural soil generally increased with an increase in applied stress. It was 
higher also far stresses applied parallel to the stratification than for stresses applied normal 
to the stratification. It is also noted that the modulus for the compacted backfill assumed values 
somewhat higher than for the natural soil. 

Variation in soil modulus for different types of tests is to he expected. The rate of applica- 
tion and duration of load undoubtedly inflmence the magnitude of the soil modulus. The consoli- 
dation, triaxial, and load-bearing tests are relatively slow tests, whereas the soniscope test 
is analogous to a rapid test wherein the load is applied for a short period of time. Other factors 
such as side wall friction for the consolidation test specimens and the distribution of stress in 
the triaxial test specimens also undoubtedly influence the test results. 

Figures 4.48 and 4.49 present typical exmples of stress-strain diagrams taken from the 
Project 3.8 report. The consolidation test results (Figure 4.49) show a characteristic slow 
take-of[ at low stress on the first compression cycle; the second cycle loading curve does ret 
display this behavior to the same deyree. This behavior might be due to some sample-apparatus 
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Figure 4.47 Normalized horizontal velocity spectra, Priscilla 
5G-foot depth, 0.5 percent critical damping. 
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Figure 4.48 Strcss-strain curves, undisturbed soil, Frenchman Flat. 
(Reprintcd from Reference 14.) 
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. increments of 10.0 psi and 2.5 psi,respec- 
Max, langent mad. of compression 914,200 ftlvely. es 2 ee aaeee 


pei. 
Max tongent mod.ol compression 219,050 


psi. 


Figure 4.49 Project 1.4 constant ratio of applied stress, triaxial tests, 
undisturbed soll, 50-foot depth. (Reprinted from Reference 14.) 
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Figure 4.50 Stress-strain diagram, Statlon 4 
(GR: 750 feet), Shot Priscilln. 


/ 
12,600 psi —>/ 


/ 


34,000 pares 


10 15 2.0 2.5 


Strain (Gage 6SVS50) , ppk 


3.0 


Figure 4,61 Stress-strain diagram, Statlon 6 (GR: 1,050 feet), Shot Priscilla, 
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coupling alignments which must be affected before the reliable measurement Is taken. Also, it 
should be noted that the tests indicate strains in excess of 4 percent for 100-psi applied stress 

in the first cyclo; the second cycle produces only about 1 percent strain for the same applied 
stress. The triaxtal results of Figure 4.49 indicate i to 1.5 percent strains for: ator stresses 
approaching 100 pst. 


4.6.2 Deduced Stress-Strain. Since stress and strain were measured independently at vari- 
ous depths in Frenchman Flat soil (Project 1.4), it would be possible to construct a measured 
stress-strain diagram. This possibility is severely lMmited by the fact that so many of the strain- 


TABLE 4.6 PROJECT 3.8 SOIL SURVEY RESULTS 


Natural Soil Compacted Backfill 
Item Ranged . Ranged 
From Ta Average From To Average 
Modulus of deformation 
{consolidation test) 
Normal to stratification 
SO psi 2,140 12,150 5,300 4,580 20,700 6,630 
100 psi 3,000 13,900 6,180 4,000 20,700 8,665 
Parallel to stratification 
50 psi 11,580 16,500 14,560 _— -- - 
100 psi 10,30 26,780 18,550 _— _ _ 
Modulus of compression 
{triaxial test) 9,050 19,050 13,300 _ — _ 
Compressive modulus 
(soniscope test) 31,000 255,000 96,200 _ _ _ 
Modulus of elasticity 
(plate bearing test) 
50 pst 6,100 6,600 6,450 _ — _ 
100 pst 4,900 3,900 4,400 — — —_ 


time records were elther lost or are in- question. However, it is possible to construct such dia- 
grams for the 50-foot depths at both Stations 4 and 6; these deduced stress-strain diagrams are 
shown in Figures 4.50 and 4.51. 

The form of the stress-strain picture of Figure 4.50 ls quite different from the Project 3.8 
curves of Figure 4.48. This is due probably, to the rather abrupt increase of stress (4CV50) 
at this station, and it results in a tangent modulus of about 87,500 psi in the early portion of the 
curve. Calculation af the modulus corresponding to the peak strain yields 17,300 psi, and the 
untoading modulus is about 32,500 psi. Also significant is the fact that the soil attains a strain 
of less than 1 percent (about 0.8 percent) at a stress of 130 psi, which means that the Project 
3.8 triaxial tests approximated the working range af dynamic soil response more exactly than 
aid the compaction tests. The residual strain (after unloading) of about 0.4 percent also agrees 
with the triaxial test residual magnitudes. 

The situation at Station 6 is also atypical but for the opposite reason from that at Station 4. 
At Station 6, rhe stress appears to increase unusually slow while the strain is attaining readable 
values. The best choice of tangent modulus corresponding to peak stress and strain computes 
to 12,600 psi, 2 value which again agrees best with the modulus obtained from the Project 3.8 
triaxial tests. The curve also indicates 2 small residual strain at this depth. 

It would be unwise to draw any firm conclusions from the small amount of data available. 
Tentatively, it can be said that, as presently conceived, the Liboratory triaxial tests are more 
useful in correlation with blast results than are the compaction tests. Still to be resolved is 
the question of the relationship between the tangent modulus so computed and the modulus deter- 
mined from measurements of seismic velocity, 
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Chopter 5 
CONCLUSIONS and RECOMMENDATIONS 


5.1 CONCLUSIONS 


5.1.1 Instrument Performance. For Operation Plumbbob Project 1.4, full length records 
were obtained on 75 percent of the total gage channels and partial records on an additional 12 
percent. The failures were caused principally by a jammed oscillograph recorder and cable 
trouble on the strain gage installations. 

Five out of the ten short-span strain gages installed on the sides of the two deep holes failed 
to give useful data, and three of the five useful records were incompicte. [It is concluded that 
the reliability of these gages could be enhanced by improving the method of coupling to the soil 
and by redesign of the cable connector at the gage end. 


5.1.2 Acceleration and Velocity Wave Forms. For ground shock accelerations induced by 
air blast, there is no ideal wave form. The integral of vertical acceleration, vertical particle 
velocity in the soll, is the more familiar of ground motion parameters since it bears a direct 
relationsh{p to the overpressure. It is possible to separate the predominant earth acceleration 
and velocity wave forms into six groups asscciated with the character of the input; namely, 

1. Air-blast wave ideal and overpressure high, traveling superseismically. 

2. Same as Item 1 above, execpt lower overpressure. 

3. Air-blast wave traveling transseismicaliy, first indicatiun of outrunning. 

4. Air-blast wave in early stages of precursor development, peak overpressures high, 
traveling super seismically. 

S. Air-blast wave extremely disturbed, in most violent precursor region, traveling super- 
or transseismically. 

6. Camplete outrunning, refracted ground wave arrives before local induced effect, air wave 
traveling subseismically. 

Even though the blast wave may be traveling with superscismic velocity, i.e., no outrunning, 
the signal produced by the refracted wave may be sufficient to modify the velocities in the later 
portions of the ground shock isturbance. At depth, outrunning occurs earlier than at the surface. 


5.1.3 Signals from Remote Sources. When the propagation velocity of the air-blast wave 
front is larger than the compression wave velocity in the soil, signals will not be observed 
underground prior to the arrival of the air wave over the gage. AS the air wave velocity de- 
creases with distance irom ground zero, the information given to the earth will eventually 
arrive before the arrival of the local signal. This may occur either when a near-surface seis- 
mic velocity is greater than the air-blast velocity or by transmission through the soil along 2 
curved path which dips into the higher-velocity lower earth strata, t.e., Seismic refraction. 

For Shot Priscilla, it was concluded that outrunning occurred at the ground surface at about 
2,500-foot ground range, duc to a signal originating at about 1,900-foot ground range. However, 
it is possible for outrunning to occur at closer ranges for deeper measurements, and refracted 
signals may be evident on records alter the arrival of the local effects. Although Project 1.4 
acceleration results (measurements to 1,350-foot ground range) showed no outrunning, they 
showed presence of refracted signals from remote sources, particularly at the deeper gages. 

The remote source disturbances, whether ovtrunning or not, are extremely hothersome when 
one ts attempting to integrate the acccleration-time records to obtain velocity-time and displace- 
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ment-time. Any criteria used for determining the time of zero velocity may be uselesy in the 
presence of remote yource signals. 


5.1.4 Attenuation of Ground Shock with Depth. Maxlmum downward acceleration data indi- 
cates that between 5- and 10-font depths the attenuation varies between 30 to 45 percent, cxcept 
at the 550- and 650-fvot ground ranges where the attenuation is practically negligible. The 
results at ercater depths, taken into consideration with the measured seismic profile at French- 
man Flat, suggest that the wave theory concerning energy transfer at an interface between two 
materials may be valld. The theory states that: 

1. if an aceclerometer is located at or near the transition from a hard to a soft material, 
the peak acceleration Ls Increased. 

2. If an accelerometer is located at or near the transition from 2 soft to a hard material, 
the peak acceleration is decreased. 

The peak outward horizontal acceleration at the 10- and 50-foot depths shows less attenuation 
with depth than does the corresponding peak downward acceleration; numerically, the attenuation 
ts about 40 percent for the outward component and about 80 percent for the downward. 

When all pertinent data are included, It is found that power law decay with depth of peak 
downward acceleration agrees best with the experimental results. 

For peak downward velocity (or velocity-jump data), an exponential decay Law appears to 
represent the data better than a power law. On Shot Priscilla, at the two stations where meas- 
urements were obtatned tn 50-foct depth, the attenuation of peak downward velocity is similar 
and is apparently little influenced by the underlying soft layer documented in the refraction 
Seismic survey. 

The peak outward Velocity data from Shot Priscilla showed somewhat less attenuation with 
depth than the downward component at 275-psi overpressure; at the larger ground range (100- 
psi overpressure), the outward peak velocity at 50-foat depth was twice that at 10 feet. This 
effect was traced to the influence of signals {rom remote sources closer to ground zero. 

The vertical velocity-time curve is similar to the overpressure-time curve when the air- 
biast-wave 1S Ssuperscismic. The velocity peak is, however, more rounded and progressively 
lags the wave front by an increasing amount as the depth increases. Below 10-foot depth, the 
rise time to peak velocity increases significantly, even for rapid air-pressure onset. 

The attenuation of peak displacement, obtained from successive integration of the acceleration- 
time data, corresponds closely to the attenuation of peak velocity. However, the larger than ex- 
pected displacements below 5 feet at the 100-psi overpressure station can be explained by a re- 
fracted ground wave Superimposed or the motion Induced by the main air-blast wave. It should = 
be noted that the double integration method of obtatning displacement-time is probably not reliable 
for times beyond the peak value. 

The attenuation of maximum vertical stress, as determined by Carlson-type gages, Is slight 
between the surface and 5-foot depth; the attenuation at shallow depths is probably a function of 
the duration of the overpressure input pulse, Lhe shorter the input pulse the greater the attenua- 
tion with depth. 

For the deeper stress measurements, the pattern of attenuation follows the rule of halving 
the peak stress for each 10-foot Increase in depth. However, the increase in stress at 50 feet 
tends to deny the validity of this rule. 

Although several possibilities are offered to explain the seeming anomalous Stress results 
at 50-foot depth, it is probable that the measurements are influenced to some degree by all the 
suggested mechanisms. That is, the stress measured at 50 feet may be a function of the soil 
layering structure, the character of the backfill, and inequalities of planting or tamping methods 
in the vicinity of the gages. 

In general, the stress-time wave forms compare well with the air overpressure input; below 
10-ffo0t depths, the time-to-peak increases noticeably over that observed on the surface air- 
blast pressure. 

It is difficult to determine the relative validity of the {rce-field stress results obtained on 
Project 1.4. The control of backfill operations during Priscilla was considered optimum, and 
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extraordinary care was taken in the placement of gages; still the stress measurements appeared 
internally inconsistent. ft Ls concluded that, unless an alternative mcasurement method can be 
devised, stress must be derived from other primary measurements such as strain and accelera- 
tion, which are more independent of hole size and character of the backfill. 

In addition to the short-span strain gage results, it is possible to compute the average verti- 
cal strain between accelerometer positions from the doubie integrations (i.e., displacement- 
time plots). At the 275-psi overpressure station, the peak strain sulfers a pronounced decrease 
between 1- and 30-foot depths. At greater depths, the vertical strain levels off to a constant 
value. On the other hand, the data obtained at the 100-psi overpressure statian Show practically 
no change in vertical strain with inerecasing depth. This difference in behavior is probably tracc- 
able to the character of the input overpressure; the rise times associated with the 100-ps! input 
are significantly longer than at the high pressure station. 


5.1.5 Ground Shock ard Overpressure. When the ratio of peak acceleration to peak overpres- 
sure is plotted versus peak overpressure, the values at high overpressurcs (above 200 psi) are 
high, decreasing as the overpressure drops and as the precursor wave develops, finally increas- 
ing somewhat below 10 psi when the precursor has dissipated. The effect of device yield, if any, 
does not appear to be discernible. No systematic variation with height of burst ts to be found; 
also, the trend toward higher ratios at extreme ground ranges may be associated with outrunning 
of the ground wave. 

In conformity with elastic theory, subsequent to outrunning of the ground wave, the velocity- 
jump peak overpressure ratios increase with decreasing pressure {i.e., increasing ground 
range). There appears to be no systematic variation with either yield, overpressure level or 
wave form, or test area. The experimental ratios agree well with the theoretical result, par- 
ticularly if the velocity of propagation of the ground wave is determined from measured accel- 
erometer arrival times. 

The fact that the finite amplitude stress waves appear to propagate sloxcr than the seismic 
waves suggests plastic flow conditions in the soil; also, relative significar™. of the soil modulus 
determined from seismic velocity is not, as yet, completely clear. 

Consideration of the peak vertical displacement-overpressure impulse ratios indicates a 
definite increase in the ratio with increasing peak overpressure in the range of 50 to 300 psi. 
However, the wide scatter in displacement data precludes any firm conclusions. Comparisons 
between the Project 1.4 peak displacement data (obtained by integration of sccelerograms) and 


the Project 1.5 long-span Strain results show good agreement between the ‘ two Independent 
measurements. 


5.1.6 Response Spectrum of Ground Motion. The response spectrum is defined as the maxi- 
mum response of a linear, single degree of freedom, spring mass system, relative to the motion 
af the ground. The displacement spectrum may be determined either by direct measurement 
with reed gages, ty numerical calculation of the response of a linear system (e.g., accelcrom- 
eter) to the ground motion, or through the use of the electrical-mechanical analogy. An accel- 
erometer-galvanometer system will yield faithful results if its natural frequency is more than 
twice the fundamental frequency observed in the ground motica. 

From the results of the displacement response Spectra, it does not appear that the change in 
character of the response can be associated with the local ground wave, but this change probably 
ts due to interference of the refracted ground-transmitted wave. 

Using overpressure positive phase duration as a normalizing factor, it is possible to show 
correlation of response spectra obtained under different input conditions. The normalized veloc- 
ity spectra for 50-foot depth when compared with spectra at 5- and 10-foot depths show similar 
maxima although the frequency at which this maximum occurs decreases with increasing depth. 

Weighing the relative advantages and disadvantayes of the reed yaiges and accelerometers fer 
determinaticn of response spectra, one concludes that neither gage result is completely suffi- 
cient for definitive spoctra; the most advantageous Inst rumentation probably will consist of both 
types placed as a pair. 
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5.1.7. Soil Stress-Strain Considerations. Since stress and strain were measured independ- 
ently (versus time) at vartouy depths in Frenchman Flat soil, it {s possible to construct 2 
measured soil stress-strain diagram. Such diagrams may be constructed for the 50-foot depths 
at Stations 4 (275-psi input) and 6 (100-psi input). 

The deduced stress-strain relations appear to ugree best with the laboratory triaxial tests 
performed upon Frenchman Flat soit samples by Waterways Experiment Station (Project 3.8). 
The consolidation test results yicld soil moduli which are much lower than those indicated by 
the deduced stress-strain curves. Tentatively, it can be concluded that the Laboratory triaxial 
tests are more uscful in correlation with blast results than are the compaction tests. However, 
Still to be resolved is the question of the relationship between the tangent modulus so computed 
and the modulus determined from mezsurements of scismic velocity. 


5.2 RECOMMENDATIONS 


It is recommended that future work on the problem of ground motion induced by air blast 
include the following: 

1. In conjunction with any future ground motion measurements, means should be provided 
fcr determination of the seismic signal arrivals at each measurement position; a geophone 
would serve this purpose. 

2. Attention should be given to devising laboratory soil test methods which give promise 
oi correlating with field test results, 

3. Work should continue for the development of a reliable gage (or gages) which will meas- 
ure directly velocity-time and/or displacement-time of ground motion. These measurements 
would almost completely eliminate the present uncertainty associated with integrated accelera- 
tion data. 

4. The theoretical aspects of the problem should be pursued vigorously, from the standpoint 
of pure elasticity and, also, to include the dissipative and layered-g6pects of the soil medium. 

5. In ground motion measurements an independent measurement of the motion of the ground 
relative to some fixed point in the earth made by a mechanical gage is useful in checking dis- 
placement results from integrations of accelerograms. 
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SYMBOLS and OVERPRESSURE WAVE FORM CLASSIFICATION 


> 


Appendix A 


SYMBOLS 


= acceleration, g (vertical) 

= density, lb/cu ft 

= Youngs" modulus, pai 

= strain, vertical 

= applied air pressure, psi 

= total rise of first peak of 
input pressure wave 

= total rise of second peak of 
input presaure wave 

= stress, vertical, pel 

= Poissons’ ratio 

= total rise time of stress (or 
velocity) wave at a particular 


depth (assuming half-cosine 
wave form of rise), 200 


= vise time of Brst peak of 
iopat pressures wave 


= rise time of second peak 
of input pressure wave 


= rise times of stress (or 
velocity) wave produced by 
travel to depth y 


= particle velocity (vertical), fps 


= seismic (compressional) wave 
velocity, fps 


= range, ft 
= depth, ft 


-y33 
@ + 9 @— 2a) 


AIR BLAST WAVE FORM, PRECURSOR TYPE I 


mon 


AIR BLAST WAVE FORM, GENERALIZED 
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A.2 OVERPRESSURE-WAVE-FORM CLASSIFICATION 


Type 


0 A sharp rise to a double-peaked maximum; 


panies aloes together tn toes and opprackaeacty peaked shock wave but is usually recorded as a 
equal in amplitude. 


Deacriptioa of Form Relation to Previous Type 
In its ideal form, it is the classical single- 


double-peaked wave. 


General Typical 


lL sA sharp rise to first low peak followed by cithe= = The first low peak indicuten the wxistence of a 


a platessa or a slight decay, then a higher second distmrbance which travels faster t= the main 
peak preceding the rapid decay. Time interval wave. This type is distinctly nonclassical. 
between first and second peaks can vary sig- 

nificantly; shock-like rises are evident. 


pees eT Men, 


Typical 


Same as Type 1 except that second peak ls less = The second peak has decayed to a lower value than 
than first. the first and has became more rounded and lass 
distinct. Second perk finally disappears. 


General Typical 
A first large. rounded maximum followed by The first peak of Type 2 has developed to become 
decay; then a later, usually amaller, second the rounded maximum, while the second peak has 
peak. Pressure rises may be slower than for decreased in magnitude with respect to the first. 


Ty 2 


Oo Ce 


Typical 
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& long-rise-time flat-topped form which ex- The relatively sharp pressure rise of Type 3 has 


hibits a lang decay Ume and much hash. deen replaced by a slew rise and the second 
peak has disappeared. 
ee 
General Typicsl 
A pressure clas to a rounded plateau which Too single-poaked hashy form of Type 4 seems 
ts followed by a slow rise to 3 second. higher to develop a compresatoa-type eecond peak, 
peak. which may be tho Gent ls dicatios of the return 


Of the resin were. 


General Tytacal 
A clear-cut double peak form with a rise to This ie clearly a cleaned-up Type 5. with the 
a piatren which elopes upward, then a shock compression-type second peaks becoming shocas. 
cise to = peak. 
op eee 
Genera! Typicat 
A shock rise to 2 peak followed by either a ‘The second peak oF Type & has overtaken tly 
slight: gentle rise, = platean, or in lus ex- first peak. reeniting in a wave form which is 
amples, a alow decay. Close to classic, sharp, single peak is ont evident. 
i ce ae Ra 
Gereral Typical 
Refers w Type T in region of regular reflec- Second rise due to reflected wave. 
tion where e second (reflected) shock front 
is eviden:. 
wR rd ~ 
8 Aciassical wave form. Starp alogie-peaked form, followed by clasaic 
decay. 
° ES (ee 
General Typical 
Refers to claseieal wave fo>= {9 regis= of Secoad rise dur to reflected wave. 


reguar reflection. 


er J Bape. 
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Appendix 8 
CATALOG of GROUND MOTION MEASUREMENTS 


TABLE %1 SUMMARY OF NEVADA SURFACE AND AIR NUCLEAR BURSTS 
ON WHICH STRONG GROUND MOTION WAS MEASURED 


Y. Yucca Flat, and FF, Frenchman Flat. 


—_et Height 
Coez=“o0n Shot Yield of Burst Date Area 
ke fe 

Buster-Jangie s 1.19 3.5 November 1951 NTS-Y 
Tumble r-Snapper 1 1.05 793 April 1952 NTS-FF 

2 2.15 1,109 April 1952 NTS-Y 

3 30 3.447 April 1952 NTS-Y 

4 19.6 1,040 May 1952 NTS-Y 

Upshot-Knothole 1 16.2 300 Marsch 1953 NTS-Y 
9 26.0 2,423 Mzy 1953 NTS-FF 
10 14.9 524 May 1953 NTS-FE 
Plumbbob Priscilla 39 700 dane 1957 NTS-FF 

Whimey 19.2 500 September 1957 NTS-Y 

Galileo 11 500 September 1957 NTS-Y 

Smokey 43 700 August 1957 NTS-¥ 

Stokes 19 1,500 Avgust 1957 NTS-Y 


Charleston 1.5 1,500 September 1957 NTS-Y 
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TABLE 8.2 NUMBER OF GROUND MOTION MEASUREMENTS, NTS 
Y, Yucca Flat; FF, Frenchman Flat; H, horizontal direction; and ¥, vertical direction. 


Operation 


dangle, ¥ 


Tumbler, FF 


Tumbler, Y 


Upshot-Knothole, FF 


Upshot-Knothola, Y 


Plumbbob, FF 


Plumbbobt, Y 


Overpressure 


psi 
>60 


>10 


Total 


>10 


Total 


>50 
>10 
Tolal 


>60 
210 


Total 


>10 


>60 


>B60 


>10 


* First source is alr blast. 


t Reed yage data. 


Surfaco 


2u 
av 


3H 
6Y 
8s 


ily 


1 


ev 


Depth, ft 

6 10 20 
TT ie 0 
ov ov 
™ il 
qv av 
Wn 4 
live o4V 

WwW 

Mv 

2 

10Vv 

WV 

wv 

6H 

30V 

2H 

8H 

2V 

13V 

2V 

MW 

ul 
4ll 

TV 12V0 OV 


18 
6v 


60 


2u 
ev 


60 


3 
dv 


aul 
av 


100 200 


iv 


Source 


of Data* 


WT-414 
WT-388 


WT-fl4 
WT- 514 


WT= 514 
WT-517 


WT= 782 
WT-T11 
WT-716 


WT-T12 
WT-716 
ITR= 1403 
ITR=1404 
ITR-1405 


ITR= b487 
ITR~ 1487 


TABLE B.3 NUMBER OF GROUND MOTION MEASUREMENTS, EPG 
H, horizontal direction; and V, vertical direction. 


Operation Depth QOvurpressure Numbor Sourco of Data® 
fe pai 
Greenhoure lto 6f >50 4H 
4vV 
>10 LiH WT-S3 
10V WT—69 
Total 1311 
13Vv 
Ivy 17 >50 1H 
iv 
>L0 SH WT —602 
3Vv WT —9002 
Tol 4H 
sv 
Castle 1S to 17 >10 Sit WT—320 
2Vv WT —9002 
Rodwing 25 > 1H ITR—1302 
6V ITR—-1364 
Bardtack 1 to 100 >” 17H ITR=—1613 
19V 


* First source is air blast. 
t Varies with station. 


TABLE 2-4 SUMMARY OF MOLE ROUNDS (256-POUND TNT) 


Height of Burst, & o 0.83 1.65 3.18 6.35 
4 (scaled height of burst for HE) 0 0.13 0.26 0.50 1.0 
Test Area 
Utah dry clay 1 1 1 2 1 
Nevada gravel and sand mix 1 I 2 2 1 
Californiz wet sand 1 1 - 
Califorata moist clay - 1 - - - 


TABLE B.S NUMBER OF GROUND MOTION 
MEASUREMENTS, JANGLE HE-4 


Height Number of 
Ctarg of Burst Careresstme Observations, 5 ft 
€ psi 
2,560-Ib/TNT 2 >50 au 
4V 
>10 8H 
8Vv 
Total 10! 
lov 
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Appendix. C 
RECORD INTEGRATION ond INSTRUMENT RESPONSE 


Ground motion instrumentaUon has been confined, 
except for the long-span displacement gages used on 
Shot Priscilla, to accelerometers. Determination of 
velocity and displacement requires numerical inte-— 
gration of these records. Prior to Shot Priscilla, 
this was usually done with a planimeter as a two-step 
process. Integration of Priscilla records has been 
performed sumerically using the IBM 650. Accelera- 
tion records were read at %4-msec {atervals through- 
out the portion of maximum fluctuation. As the records 
smoothed out, this time interval was lengthened suc- 
cessively to 1, 2, and 4 msec, mainly as an economy 
measure. Records were read to the end of the alr- 
blast.positive phase or past the point of apparent sig- 
nal, whichever was longer. 


C.1 ACCELERATION SASELINE SHIFT 


It Ls an indisputable fact that the ground motion due 
to air-blest pressure must cease at some time alter 
the passage of the blast wave. Most integrations indi- 
cate, in varying degree, that the velocity at the end 
of the integration is not zero. This result can be 
interpreted as an acceleration record baseline shift, 
which is affected during recording or during reading 
of the record. If it can be assumed that there are no 
frequency response problems (sea Section C.2), the 
source of the difficulty can be traced to the character 
of the acceleration-time wave form. Reference to 
Figure 4.12 shows tkat the duration of the first accel- 
eration peak Is a small fraction of the total record 
length. Therefore, a small error, perhaps only one 
or two percent of the peak acceleration, will accumu- 
late 23 time increases to result in a significant error 
in velocity at the end of the integration. 

Speaking instrumentally, there are many possible 
ways the acceleration record baseline (zero signal 
reference) could shif. However, in this analysis the 
shift is considered as a reading error, a conclusion 
which Ls substantiated by the fact that the amount of 
correction that is necessary to achieve zero velocity 
at the end of the record is frequently within the least 
count of the reading equipment (see Table C.1}- 

Figure C.1 shows some typical exampics of velocity- 
time wave forms obtained from types a, b, or ¢ accel- 
eration-time wave forms (Figure 4.1). It is thought 
that a good eriterton for zero velocit; in such cases 
is that the velocity equals zero at a time after signal 
arrival equal to the duration of the air-blast inpuc 
positive phase. For this criterion, Figure C.1 indi- 
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cates how the baseline would be shifted on the velocity- 
time plots; subsequently, the velocity-time plots rela- 
tive to the shilted baselines would be integrated to 
obtain the corrected displacement-time plots. It should 
be pointed out that, in the case of gages buried deeper 
than SO feet in a soil whose seismic velocity varies 
significantly with depth, it would be unwise to apply 

the above criterion for zero velocity even for Ideal 
overpressure inputs. Each record would have to be 
regarded as a special case. 

As the peak acceleration decreases, the problems 
of accurate record reading diminish with the exception 
of greund waves trom remo sources, whith cesult in 
exceptionally complicated records. A typical example 
(type f) is shown in Figure C.2, where the upper por- 
tion of the figure illustrates the original integration 
without correction. Beyond 800 mscc the velocity 
increases (negatively) linearly. The iinear baseline 
shift required to correct the velocity to zero at the 
end of the greind motion is shown as well as the cor- 
rected velocity. However, it is at once apparent that 
ths criterion of zero velocity at the ond of the air-blast 
pressure pulse is not applicable when signals from 
remote sources are present. Far these, each uncor- 
rected velocity must be plotted and a best judgment 
made as to the time zero velocity is attalned; thus It 
is Important that the acceleration-time record be read 
to times beyord which the signal apparently has settled 
down to small amplitudes. In complex cases, the prac- 
tice has been adopted that the baseline be adjusted to 
maximize the resuiting peak displacement. The ais- 
placement-time records will then indicate a residual 
displacement equal to the peak displacement. If more 
than one choice of baseline shift appears reasonable, 
the extremes are calculated. 

Uf the velecity records appear to require other than 
a single linear correction; e.g-, a series of Linear 
corrections attached end-to-end, the tata are consid- 
‘red suspect and not corrected or tabulated. 

Table C.1 includes a summary of the magnitude of 
the corrections made in the integrations of the Pris- 
¢cilla aeceleration-time records. In terms of the peak 
acceleration, the baseline shifts varied between 0.06 
and 4.4 percent, which corresponds to between 0.1 and 
5 counts on the electromechanical reader (Oscar). 

The reader is set up so that one count equals slightly 
less than Yygo-inch deflection on the gage record and 
the reproducing accuracy of a single reader is approxi- 
mately + 2 counts. Thus it apoears that the need for a 
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TABLE C.1 BASELINE CORRECTIONS, SHOT PRISCILLA 


(1) and (2) indicate alternata baseline choices. 


o av 4D 
Gage AUncor A Cor K Gor co VUncor V Cor Vv Gor DUncor DCor D Cor 
4 g pct sec sec pet t ft pet 
1v5 363.7 387.7 1.68 a 26.7 25.3 314 214 0.693 208. 
LV1O 22.9 225.9 1.78 3 13.9 14.8 Gil 0.298 0.343 13.1 
2V5 178.9 173.4 3.17 4.5 20.7 17.9 14.98 1.47 0.586 180.2 
2vV10 150.6 151.9 0.87 L5 19.6 20.2 235 0.718 0.964 25.5 
3V5 120.5 120.9 0.27 6.3 16.1 16.3 Lis 0.717 0.839 14.5 
3V10 123.6 129.3 0.22 1 13.3 13.1 1s 0.831 0.555 13.7 
4V1 189.3 188.2 0.56 1 17.0 15.9 7.20 0.898 0.664 33.2 
4V5 $3.9 $4.7 1.35 2 14.4 18.3 &53 0.543 0.768 23.3 
4vV10 38.5 39.2 1. 2 10.1 10.9 7.13 0.373 0.521 28.4 
4V20.1 4180 6.6 1.35 6 6.91 7.64 9.65 0.324 0.486 3295 
4V30A 37.39 38.1 0.35 2 7.72 7.92 254 0.345 0.2388 1L1 
4V50 13.2 23.5 2.27 3 4.51 4.96 3.13 0.221 0.311 28.9 
4VS0A 125 12.8 2.16 5s 4.23 4.66 8.73 0.207 0.2489 28.4 
4H10 14.2 14.3 1.06 2 1.27 1.42 3L1 0.016 0.073 73,5 
4850 5.08 $.07 0.10 0.3 0.880 0.885 0.56 0.014 0.014 —_ 
sv5 24.3 _ _ _ _ _ _ _ _— —_ 
svV10 16.7 _ — —_— _— _ _ — — _ 
$¥1 (2) 21.8 21.5 L233 1 6.72 6.12 9.87 0.300 0.218 37.6 
(2) 2L3 21.7 0.34 0.1 6.72 655 2.83 0.348 0.299 16.4 
6V5 (1) 16.7 16.7 0.08 0.1 6.05 6.02 0.53 0.274 0.267 262 
(QB 16.7 16.8 0.33 0.4 6.05 6.18 217 0.273 0.307 1L1 
6V20 (1) 10.5 10.5 0.23 0.2 4.13 4.07 Ls5 0.322 0.306 §.23 
(2) 10.5 10.5 0.20 0.2 4.43 4.08 1.37 0.323 0.307 5.21 
6V20A (1) 9.80 9.93 L.32 2 3.6% 3.96 B.76 0.221 0.275 19.6 
{23 9.80 9.97 1.68 z 3.61 4.06 10.9 0.220 0.293 24.9 
6V20 6.49 6.53 0.54 0.5 3.49 3.59 278 0.191 0.214 10.8 
SV30A 6.32 6.40 L2e8 3 3.32 3.56 6.67 0.166 0.209 20.6 
6V50 (1) 262 2.72 3.71 4 L397 2.30 14.2 0.125 0.169 26.0 
(2) 2.62 274 4.38 4 1.97 2.36 16.4 0.226 0.182 30.8 
6H10 1.30 1.31 0.53 0.3 0.383 0.400 4.25 0.006 0.008 25.0 
6HSO 253 254 0.71 2 0.877 0.920 4.67 0.020 0.022 9.09 
75 9.23 9.16 0.85 1 3.09 2.84 8.55 0.376 0.239 57.3 
TV10 (1) 4.84 484 0.06 0.2 1.50 1.49 0.80 0.086 0.090 4.44 
(2) 4.84 4.78 1.22 2 1.50 1.67 20.11 0.178 0.154 15.6 
Average 1.17 Ls 6.25 32.20 
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Figure C.1 Typical velocity-time wave forms from Types 2, 3, ore 
acceleration-time wave forms. 
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Figuee C.2 Uncorrected and corrected vertical velocities. 
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baseline shife in the Integration process can be traced, 
with few exceptions, to reading errors. This fact puta 
a great emphasis upon accurate record reading, if 
point-by-point integration Ls to be attempted. It Is 
also apparent from the table that the corrections In 
terms of average percent of peak amplitude were 
smallest for acceleration (1.17 percent), larger for 
velocity (6.25 percent), aod largest for displacement 
(32.2 percent). 

in some cases, particularly on records obtained at 
Stations 6 and 7, two independent baseline corrections 


PEAK NEGATIVE ACCELERATION-9g 


100 


Figure C.3 Comparison of ERA 


wave, for damping between 0.5 and 0.6 eritical, tran- 
sient pulses are reproduced [alrly well for pulye fre- 
quencies of one sixth the gage frequency. When the 
pulse frequency becomes equal to or greater than one 
half the gage frequency, considerable error in gage 
response along with phase distortion may be expected. 
For the half sine pulse, the errors are approximately 
20 percent. If the damping ts only 0.4 critical, over- 
shoots of 50 to 100 percent occur. 

In general, the frequency response of accelcrom- 
eters used in weapons cffects experiments has in- 


1000 
RANGE -ft 


and Wiancko accelerometer data, 


Tumbler Shot 1, $-foot depth (numbers adjacent to data indicate 
accelerometer natural frequencies). 


were made and the integrations carried througn using 
both. Deviations in peak velocity from the two choices 
varied between 2.5 and 10.8 percent, whereas the peak 
displacements were separated by as much as 41.5 per- 
cent. The largest deviations between the two bascline 
choices were observed on the 7V10 record; it is prob- 
ably more than coincidence that this record displays 
the most complex wave form due lo ground shock sig- 
nals from remote sources. This result leads to the 
conclusion that integration of acceleration-time records 
possessing complex wave forms invo!ves a good deal of 
judgment and is necessarily less accurate than integra- 
tion of simple records, f.e., those with local effects 
only. 


C.2 INSTRUMENT RESPONSE 


The complex pattern of the ground accelerations 
makes a precise statement of error due to frequency 
response intractable. A few general remarks can be 
made, however, for simple inputs such as a half siac 


creased markedly since Operation Greenhouse. The 
natural frequency of accelerometers used on Operation 
Jangle varied between 10 and 140 cs, (Reference 19), 
on Operation Tumbler between 30 and 196 eps, (Refer- 
ence 2), 2nd on Operation Upshot-Kne:bole all vertical 
accelerometers had frequencies of 450 cps, (Referencs 
9). In the latter, the frequency response was limite? 
by the galvanometer circuit at 2 slightly lower vaiue. 
Shot Priscilla response characteristics were similar 
to those af Operation Upsbot-Knaothole. 

A graphic illustration of inadequate frequency re- 
sponse is found in the comparison of ERA accelerom- 
eter results (Reference 28) and Wiancko accelerometer 
tesults on Tumbler (Reference 2). ERA accelerometers 
with natural {rcquencles of the order of 40 to 50 cps 
were used to back up the orimary Instrument line. 
Figure C.3 compares the two sets of data on Tumbler 
Shot 1. For accelerations greater than 5 g, the ERA 
equipment [fails to give satisfactory results. 
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Appendix O 


SCALING of POSITIVE OVERPRESSURE 


IMPULSE and DURATION 


Figures D.1 and D.2 present the major portion of 
available data on impulse and duration for overpres- 
sures greater than 10 psi-e As a guide in extrapolating, 
previously presented correlations, (Reference 29), 
results of the theoretical solution for the point source 
explosion in real air (Reference 30) are included in the 
figures, modified by 2wW theory (surface bursts). 

Impulse data on Prisctlla are in agreement with the 
theoretical curve at high pressures. At intermediate 
pressures (109 to 100 psi) measured Impulae ia from 0 
to SO percent greater (duc to precursor formation) than 
2W theory predicts. While it is understood that 2W 
theory does not apply theoretically to air bursts, it is 
believed thar impulse resulting from 2W theory will be 
a lower limit. Impulse predictions used throughout 
the report aru based on this postulate, applying a SO 
percent correction as an upper limit. 

Data on duration do nol follow the theoretical calcu- 
lations for overpressures greater shan 100 psi. It is 
suspected that this may be a limitation of the Instru- 
mentation. In order that the maximum overpressure 
does not over-range the gaya or recording instrunren- 
tation, the sensitivity of the system must be reduced 
at high pressures. A check of the SRI system indicates 
that when consideration is given to the absolute accuracy 
of reading camera records, the accuracy cf reading 
pressure amounts to approximately 5 percent af the 
peak overpressure. Using this criterion, the theoret- 
ical duration curve has been constructed for the dura- 
tion at which the overpressure reached 5 percent of its 
peak value. This curve deviates markedly from the 
theoretical phase duration (associated with zero over- 
pressure) 23 the pressure increases. However, over 
95 percent of the total positive overpressure impulse 
ia included before the overpressure drops to 5 percent 
af the peak overpressure. Hence, the theoretical curve 
of Figure D.2 is belleved to be a legitimate guide in the 
high pressure region provided alluwance 1s made for 
the increased duration of the positive phase in the 10- 
to 100-psi region due to precursor formation. 


131 


me SE A OER SPY RT ANS RR OEE pS EOL CLE OT EE 
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Figure D.1 Positive overpressure impulse versus peak overpressure. 
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